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Abstract 

Despite numerous mass extinctions in the Phanerozoic eon, the overall trend in 
biodiversity evolution was not blocked and the life has never been wiped out. Almost 
all possible catastrophic events (large igneous province, asteroid impact, climate change, 
regression and transgression, anoxia, acidification, sudden release of methane clathrate, 
multi-cause etc.) have been proposed to explain the mass extinctions. However, we should, 
above all, clarify at what timescale and at what possible levels should we explain the mass 
extinction? Even though the mass extinctions occurred at short-timescale and at the species 
level, we reveal that their cause should be explained in a broader context at tectonic timescale 
and at both the molecular level and the species level. The main result in this paper is that 
the Phanerozoic biodiversity evolution has been explained by reconstructing the Sepkoski 
curve based on climatic, eustatic and genomic data. Consequently, we point out that the P-Tr 
extinction was caused by the tectonically originated climate instability. We also clarify that 
the overall trend of biodiversification originated from the underlying genome size evolution, 
and that the fluctuation of biodiversity originated from the interactions among the earth's 
spheres. The evolution at molecular level had played a significant role for the survival of life 
from environmental disasters. 
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RESULTS 



Let us go back to the early history of our planet, and gaze at these just originated lives. They 
seemed so delicate, however they were indeed persistent and dauntless. They had a lofty aspiration 
to live on until the end of the earth; otherwise the rare opportunity of this habitable planet in 
the wildness of space may be wasted. Their story continued and was recorded in the big book 
of stratum. This story was so magnificent that we were moved to tears time and again. Was 
the life just lucky to survive from all the disasters, or innately able to contend with any possible 
challenges in the environment? Before answering this question, we should explain the evolution 
of biodiversity by appropriate driving forces. 

Again, let us go back to mid nineteenth century, and size up the situations for the founders of 
evolutionism. They were completely unaware of the molecular evolution; they knew little about 
the marine regression or transgression and paleoclimate; and they possessed poor fossil records. 
However, they still pointed out the right direction to understand the evolution of life by their keen 
insight. What is the mission then for contemporary evolutionists in floods of genomic and stratum 
data? Can we go a little further than endless debates? 

The Sepkoski curve based on fossil records indicates the Phanerozoic biodiversity evolution [H] 
ED ED, where we can observe five mass extinctions, the background extinction, and its increasing 
overall trend. The main purpose of this paper is to explain the Sepkoski curve by a tectono-genomic 
curve based on climatic, eustatic (sea level) and genomic data. We propose a split scenario to study 
the biodiversity evolution at the species level and at the molecular level separately. We construct a 
tectonic curve based on climatic and eustatic data to explain the fluctuations in the Sepkoski curve. 
And we also construct a genomic curve based on genomic data to explain the overall trend of the 
Sepkoski curve. Thus, we obtain a tectono-genomic curve by synthesizing the tectonic curve and 
the genomic curve, which agrees with the Sepkoski curve not only in overall trend but also in 
detailed fluctuations (Fig 1): 

CurveS epko ski « Curve JTectonoGenomic. 

We observe that both the tectono-genomic curve and the Sepkoski curve decline at each time of 
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the five mass extinctions (O-S, F-F, P-Tr, Tr-J and K-Pg). The growth rates of the tectono-genomic 
curve and the Sepkoski curve also coincide with each other. Hence, we show that the biodiversity 
evolution is driven by both the tectonic movement and the genome size evolution. The main steps 
in constructing the tectono-genomic curve are as follows. 

(1) We obtained the consensus climate curve (Curve _CC), the consensus sea level curve 
(Curve JS L) and the biodiversification curve (Curve -BD) to describe the Phanerozoic climate 
change, sea level fluctuation and biodiversity variation respectively (Fig 2a). (i) We obtained 
Curve JCC by synthesizing the following three independent results on Phanerozoic climate change 
in a pragmatic approach (Fig Sla): Berner's atmosphere CO2 curve 0), the Phanerozoic global 
climatic gradients revealed by climatically sensitive sediments B21 O, and the Phanerozoic 
87 5r/ 86 5r curve 0; (ii) We obtained Curve JSL by synthesizing the result in ref. [HI and the 
results in ref. [|9) [fTOl (Fig. Sic); and (iii) We obtained Curve JiD based on fossil record (Fig. 2d). 

(2) We calculated the correlation coefficients r£ v among Curve _CC, Curve JS L and Curve JiD 
(Table 1). The correlation coefficient between Curve -BD and Curve JSL in the Phanerozoic eon 
is r™ c = 0.564, which generally indicates a same phase between Curve -BD and Curve JSL. The 
correlation coefficients between Curve -BD and Curve JCC, and between Curve JSL and Curve JCC 
in the Paleozoic era are r p BC = 0.1 14 > and r p cs = 0.494 > respectively, which generally indicate 
the same variation pattern (or the same phase) of Curve JCC with Curve JiD and Curve JSL in the 
Paleozoic era. While the correlation coefficients between CurveJBD and Curve JCC, and between 
Curve JSL and Curve JCC in the Mesozoic era are r% c = -0.431 < and r^ s = -0.617 < 
respectively, which indicate a "climate phase reverse event" from same phase to opposite phase 
in P-Tr boundary. In the supplementary methods, we confirm the reality of such a "climate phase 
reverse event" by verifications for 10 group curves based on candidate climate, biodiversity and 
sea level data. Therefore, when constructing the tectonic curve based on Curve JS L and Curve JCC, 
we chose a positive sign for Curve JSL throughout the Phanerozoic eon; and we chose a positive 
sign for Curve JCC only in the Paleozoic era, but a negative sign for Curve JCC in the Mesozoic 
and Cenozoic eras (Fig Sle). 
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(3) The overall trend in biodiversity evolution is about an exponential function IfTTTl : N genus = 
Ngenus ex P(~t/ T BD)- Based on the relationship between certain average genome sizes in taxa 
and their origin time, we found that the overall trend in genome size evolution is also an 
exponential function [[[2]| [JT3J1 (Fig 3a): N genome = N° genome exp(-f/r GS ). The log-normal genome 
size distributions (Fig S2a, 3b) and the exponential asymptotes of the accumulation origination 
and extinction number of genera (Fig 2d) also indicate the exponential growth trend in genome 
size evolution. We found that the "e-folding" time of the biodiversity evolution t bd = 259.08 
Million years (Myr) is approximately equal to the "e-folding" time of the genome size evolution 
t gs = 256.56 Myr (Fig 3d): 

Hence, we can explain the overall trend in biodiversity evolution by constructing the genomic 
curve based on t gs . 

In the split scenario, we can explain the declining Phanerozoic background extinction rates 
[fT4l [fT5l according to the equation: 

rate 0+e = exp(-k GS ■ (-t + 542.0)) • rate-essential, 

where the declining factor exp(-k GS ■ (—t + 542.0)) is due to the increasing overall trend in genome 
size evolution (Fig 2c). The underlying genomic contribution to the biodiversity evolution prevents 
the life from being completely wiped out by uncertain disasters. 

So far, we have explained the declining background extinction rates and the increasing overall 
trend of the Sepkoski curve. The remaining problem is to explain the mass extinctions. Since we 
have successfully fulfilled the tectono-genomic curve to explain the Sepkoski curve, the reasons 
that caused the fluctuations in the tectono-genomic curve are just what caused the mass extinctions. 
We should emphasize here that the fluctuations in the tectono-genomic curve have nothing to do 
with the fossil data. According to the methods in constructing the tectono-genomic curve, we 
conclude that the mass extinctions were caused by both the sea level fluctuations and the climate 
changes. We refer it as the tectonic cause of the mass extinctions, which rules out any celestial 
explanations. 
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Furthermore, we point out that the greatest P-Tr extinction uniquely involved the climate 
phase reverse event, which occurred not just coincidentally with the formation of Pangaea and the 
atmosphere composition variation 0) [|T6l [TP71 . The fossil record indicates a two-stage pattern at 
the Guadalupian-Lopingian boundary (GLB) lfT8ll [fT9l Il2~0ll and at the Permian-Triassic Boundary 
(PTB) £T) 031. In detail, it also indicates a multi-episode pattern in the PTB stage E3 IE31 . 
The P-Tr mass extinction was by no means just one single event. The multi-stage/episode pattern 
can hardly be explained by the large igneous province event [|25l ||26l . We can explain the above 
two stages by two sharp peaks observed in dJCC (the variation rate curve of Curve _CC) at GLB 
and PTB respectively, which show that the temperature increased extremely rapidly at GLB and 
decreased extremely rapidly at PTB (Fig 2b). The different climate at GLB and at PTB resulted in 
different extinction time for Fusulinina (at GLB) and Endothyrina (at PTB). 

At last, we will focus on the genomic contribution to the biodiversity evolution. We can obtain 
both the phylogenetic tree of species (Fig S3a, 4c by M ci ) and the evolutionary tree of 64 codons 
(Fig 4a, S3b by M codon ) based on the same codon interval correlation matrix A. This is a direct 
evidence to show the close relationship between the molecular evolution and the biodiversity 
evolution. On one hand, the result is reasonable in obtaining the tree of species. This universal 
phylogenetic method based on M ci applies for Bacteria, Archaea, Eukarya and virus. On the 
other hand, the result is valid in understanding the genetic code evolution G71l Il28ll [|29l . And 
an average codon distance curve Barrier based on M codon reveals a midway "barrier" in the 
genetic code evolution (Fig 4b, S3c). Moreover, we can testify the three-stage pattern (Basal 
metazoa, Protostomia and Deuterostomia) in Metazoan origination QUI according to the genome 
size evolution. Favorable phylogenetic trees can also be obtained by the correlation matrices M gs 
based on genome size data (Fig 3c, S2c, S2d). 
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METHODS 



1 Data resources and notations 

1.1 Data resources 

(1) Phanerozoic climate change data: ref. [4], [5], [6], [7]; 

(2) Phanerozoic sea level fluctuation data: ref. [8], [9], [10]; 

(3) Phanerozoic biodiversity variation based on fossil records: ref. [1], [2], [3]; 

(4) Genome size databases: Animal Genome Size Database lf3~TTh Plant DNA C-values Database 
123; 

(5) Whole genome database: GenBank. 

1.2 Notations 

Sepkoski curve 
teconto-genomic curve 
time 

biodiversity curves 
sea level curves 
climate curves 
correlation coefficients 
climate phases 
genome sizes 
biodiversity variation rates 
derivative curves 
overall trends 
e-folding times and growth rates 
genomes size distributions and matrices 
codon interval distributions and matrices 
genetic code evolutionary curves 



Curve __S epkoski 
Curve _T ectonoGenomic 
t,T 

Curve-BD, BD, Total-BD 
CurveS L, S 1 , S 2 , S w 
Curve -CC, C ,C,C, C w \, C W 2, C^ 
r^ v ,R + ,R',AR, Q, Q',AQ 

cpi, cpii, cpiii 

G, G S p, G mean _i g, G sc [_i g, G 
ratejori, ratejzxt, rate -essential 
dJCC, dSL, dJ$D 
OT-BD, OT-GS 
tbd, kso, tgs j kcs 

M c i, M CO( j on 

Barrier, Hurdle. 
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1.3 Math notations 



Let sum(V), mean(V), std(V), log(V) and exp(V) denote respectively the summation, mean, stand 
deviation, logarithm and exponent of a vector V(i), i = 1,2, i m : 

sum(V) = J]v(i) (2) 
meaniV) = — sum{V) (3) 



std(V) = ^mean((V - mean(V)) 2 ) (4) 

logiy) = [io ge (y(i)),iog e (y(2)),...,iog e (y(/ m ))] (5) 

expiV) = [exp(y(l)),exp(y(2)),...,exp(yO ffl ))]. (6) 

Especially, let nondimiV) denote the operation of nondimensionalization for a dimensional 
vector V, 

nondim(V) = (V - mean(V))/std(V). (7) 

In this paper, we obtain respectively the dimensionless vectors Curve -BD, Curve _CC, Curve JSL, 
etc. after nondimensionalization based on the dimensional raw data of biodiversity curve, climate 
curve and sea level curve in the Phanerozoic eon. 

Let corrcoef(V, U), max(V, U), min{V, U) and [V, U] denote respectively the correlation 
coefficient, maximum and minimum of a pair of vectors V(i) and U{i) (i =1,2, i m ): 

- mean(V))(U(i) - mean(U)) 
corrcoef(V,U) = 1=1 (8) 

yjlJ^iVU) - mean(V)f Jl^m) - mean{U)f 
max(V, U) = [max(V(l), t/(l)), max(V(2), t/(2)), max(V(i m ), U(i m ))] (9) 

min(V,U) = [min(V(\),U(\)),min(V(2),U(2)),...,min(V(i m ),U(i m ))]. (10) 

Let f t (V) denote the discrete derivative of V(t) with respect to time t: 

d _ dV dV 

j t (V)-[—\ t=til) ,...,—\ t=t(i J, (11) 
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where V(t) = [V(l),V(2), ...,V(i m )] is an z m -element discrete function of time t = 
1/(1), t(2), t(i m y\. The linear interpolation of V is denoted by: 

[V(l), V(2), V(?J] = interp([t(l), t(i m )l [V(l), V(i m )l [t(l), t(i'J]). (12) 

The concatenation of function V(t) between period t([Pi]) = [t(ii),t(i\ + 1), ...,t(i 2 )] and period 
t([P2\) = [t(t2 + 1), t(i 2 + 2), t(i 3 )] is denoted by: 

[V([Pi]), V([P 2 ])] = [V(h\ V(i 2 ), V(i 2 + 1), V(i 3 )], (13) 

where Pi = [h,h + 1, i 2 ] and P2 = [i 2 + l,i 2 + 2, z 3 ] are parts of the indices. For a i,„-by-j m 
array M(i, j), let M(z, :) denote 

M(i, :) = [M(/, 1), M(i, 2), M(/, j m )]. (14) 

2 Understanding the Sepkoski curve through the 
tectono-genomic curve 

The Phanerozoic biodiversity curve has been explained in this paper. We propose a split scenario 
for the biodiversity evolution: 

Biodiversity evolution = Tectonic contribution + Genomic contribution. (15) 

We construct a tectono-genomic curve based on climatic, eustatic (sea level) and genomic data, 
which agrees with the Phanerozoic biodiversity curve based on fossil records very well. We 
explain the P-Tr extinction by a climate phase reverse event. And we point out that the biodiversity 
evolution was driven independently at the species level as well as at the molecular level. 

3 The overall trend of biodiversity evolution 
3.1 Motivation 



A split scenario is propose to separate the Phanerozoic biodiversity evolution curve into its 
exponential growth part and its variation part. 
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3.2 The exponential outline of the Sepkoski curve 



The Phanerozoic biodiversity curve (namely the Sepkoski curve) can be obtained based on fossil 
records. We denote the Phanerozoic genus number biodiversity curve in ref. [2] after linear 
interpolation by (Fig 1): 

CurveS epkoskiit) : ref. [2], (16) 

which is a 5421 -element function of time t, from 542 million years ago (Ma) to Ma in step of 0.1 
million of years (Myr): 

t = MD,f(2),r(3),...,r(5419),f(5420),r(5421)] 

= [542.0, 541.9, 541.8, ...,0.2, 0.1,0]. 1 } 

The outline of Curve _S epkoski(t) is an exponential function: 

Ngenusit) = N° genus eXp(-?/T BD ), (18) 

where the genera number constant is N genus = 2690 genera, and the "e-folding time" of the 
biodiversity evolution is r BD = 259.08 Myr. 



3.3 The split scenario of the Sepkoski curve 



We define the total biodiversity curve Total-BD in the Phanerozoic eon by the logarithm of 
Curve _S epko ski: 

Total-BD = log(Curve_S epkoski(t)), (19) 

which is also a 5421 -element function of time t. According to the linear regression analysis, the 
regression line of Total-BD on t is defined as the overall trend of total biodiversity curve: 

OT-BD = log(N genus (t)) r2m 
= k BD -(-t) + log(N° gems ), { ) 

where the growth rate of biodiversity evolution, namely the slope of this regression line, is k BD = 

1/t bd = 0.0038598 Myr 1 . 

We propose a "split scenario" in observing the Phanerozoic biodiversity evolution by separating 
the Sepkoski curve into its exponential growth part and its variation part. In this scenario, the total 
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biodiversity curve Total-BD can be written as the summation of its linear part OT-BD and its net 
variation part BD (Fig. 2d): 

Total-BD = OT-BD + BD. (21) 



Hence, we obtain the biodiversity curve Curve JD after nondimensionalization of BD: 

Curve _BD = nondim(BD). (22) 

4 The tectonic cause of mass extinctions 

4.1 Motivation 

We construct the tectonic curve based on the climatic and eustatic data in consideration of the 
phase relationships among Curve -BD, Curve -CC and Curve JSL. 

4.2 The consensus climate curve 

We denote the three independent results on Phanerozoic global climate in ref. [5] [6], [7], [4] as 
Cq, Cq, Cq respectively after linear interpolation: 

C l (t) : ref. [5] [6], (23) 

C 2 (t) : ref. [7], (24) 

C 3 (t): ref. [4]. (25) 

The missing 87 S r/ 86 S r in ref. [7] in lower Cambrian are obtained from ref. [031 for Cq. We obtain 
three dimensionless global climate curves after nondimensionalization: 

C\t) = nondim(C l (t)), (26) 

C 2 (0 = nondimiClit)), (27) 
C\t) = nondim(C\{t)). (28) 
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Hence, we obtain the consensus climate curve Curve JCC by synthesizing the above three 
results C\ C 2 and C 3 (Fig. SI a): 

Curve JCC = nondim((C l + C 2 + C 3 )/3). (29) 



4.3 The consensus sea level curve 

We denote the Phanerozoic sea level curves in ref. [8] and in ref. [9] [10] as 5 ( , and 5^ (via linear 
interpolation) respectively: 

S l (t) : ref. [8], (30) 
S 2 (t): ref. [9] [10]. (31) 
And we obtain the dimensionless sea level curves after nondimensionalization: 

S l (t) = nondim(S l (t)), (32) 
S 2 (t) = nondim(S 2 (t)). (33) 

Hence we obtain the consensus sea level curve Curve SL by synthesizing the two results S 1 
and S 2 (Fig. Sic): 

Curve JS L = nondim((S 1 + S 2 )/2). (34) 

We can obtain the derivative curves dJCC, dJ>L and d-BD respectively as follows (Fig. 2b): 

dJCC = %{Curve_CC) (35) 
at 

dSL = ^{Curve_SL) (36) 
at 

d 

d_BD = —(Curve-BD). (37) 
at 

4.4 Correlation coefficients among Curve JCC, Curve JS L and Curve ^BD 

So far, we have obtained the first group (n = 1) of curves Curve JCC, Curve JS L and Curve -BD to 
describe the Phanerozoic climate, sea level and biodiversity. They are all 5421 -element functions 
of time t. 
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There are three eras (Paleozoic, Mesozoic and Cenozoic) in the Phanerozoic eon, the time t in 
the Phanerozoic eon can be concatenated as follow: 

t=[t([P]),t([M]),t([C])l (38) 

where the indices for the Paleozoic, Mesozoic and Cenozoic are as follows respectively: 

P = [(5421 -5420), ...,(5421 -2510)], for Paleozoic from 542.0 Ma to 251.0 Ma, (39) 
M = [(5421 -2510+ 1), (5421 -655)], for Mesozoic from 251.0 Ma to 65.5 Ma, (40) 
C = [(5421 - 655 + 1), 5421], for Cenozoic from 65.5 Ma to today. (41) 

Similarly, we define the indices for the other periods as follows: 



PMC : 


for Phanerozoic from 542.0 Ma to Ma, 


(42) 


PM : 


for Paleozoic and Mesozoic from 542.0 Ma to 65.5 Ma, 


(43) 


MC : 


for Mesozoic and Cenozoic from 251.0 Ma to Ma, 


(44) 


P\L : 


for Paleozoic except for Lopingian from 542.0 Ma to 260.4 Ma, 


(45) 


L : 


for Lopingian from 260.4 Ma to 251.0 Ma, 


(46) 


L.M.Tr : 


for Lower and Middle Triassic from 251.0 Ma to 228.7 Ma, 


(47) 


M\L.M.Tr : 


for Mesozoic except for Lower and Middle Triassic from 228.7 Ma to 65.5 Ma 


■ (48) 



We can calculate the correlation coefficients r^ v among Curve -CC, Curve _S L and Curve -BD 
in certain periods respectively (Data_2): 

r^ v = corrcoe /(curve fi([p]), curve v([p])) (49) 

where the subscripts 

H,v = C,S,B (50) 
for the curves Curve -CC, Curve JS L and Curve -BD respectively, and the superscript 

p = P,M, C, PMC, PM, MC, P\L, L, L.M.Tr, M\L.M.Tr (51) 

for the corresponding periods respectively. 
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Note: The correlation coefficients generally agree with one other in the calculations between 
Curve J3D and any of Curve JSL, S 1 , S 2 , or between Curve J3D and any of Curve JCC, C 1 , C 2 , C 3 , 
i.e. in general: 

r° MV (n) ~ <vOA n,n' = 1,2, ...,10. (52) 
Therefore, the phase relationship of Curve _CC, Curve JSL and Curve J$D is generally irrelevant 
with the weights in obtaining Curve JCC and Curve JSL. The correlation coefficients are also 

irrelevant whether we nondimensionalize the curves, for instance: 

corrcoef((S '([/>]) + S 2 ([P]))/2, BD([P])) 
= corrcoef(nondim((S \[P]) + S 2 ([P]))/2), nondim(BD([P]))) 
= corrcoe f (Curve JSL([P]),Curve_BD([P])) ( } 

= r p 

SB' 

Note: The first group (n = 1) of curves Curve JCC, Curve _SL and Curve _BD is the best among 
the 10 similar groups of curves to describe the Phanerozoic climate, sea level and biodiversity. 

4.5 Three climate phases 

We propose three climate patterns CP I, CP II and CP III in the Phanerozoic eon based on the 
positive or negative correlations among Curve JCC, Curve JS L and Curve J3D. Interestingly, the 
time between the positive correlation periods and the negative correlation periods agree with the 
Paleozoic-Mesozoic boundary and the Mesozoic-Cenozoic boundary. 

(1) We have 

r p SB = 0.5929 >0 (54) 

r p BC = 0.1136 >0 (55) 

r p s = 0.4942 >0 (56) 

which indicate the positive correlations among Curve JCC, Curve JS L and Curve J3D in the 
Paleozoic era. This is called the first climate pattern (CP I); 



(2) We have 



rf B = 0.9054 >0 (57) 



13 



r% c = -0.4308 <0 (58) 
r™ s = -0.6171 < (59) 

which indicate the negative correlations between Curve JCC and Curve SL and between Curve -CC 
and Curve -BD, and the positive correlation between Curve JS L and Curve -BD in the Mesozoic era. 
This is called the second climate pattern (CP II); 

(3) We have 



r 



C SB = -0.8314 <0 (60) 



C BC = -0.8814 <0 (61) 



r 

r c cs = 0.9501 > (62) 

which indicate the negative correlations between Curve JCC and Curve -BD and between Curve SL 
and Curve -BD, and the positive correlation between Curve _S Land Curve JCC in the Cenozoic era. 
This is called the third climate pattern (CP III). 

We define the average correlation coefficient R + in the positive correlation periods: 

P P P P P P MM C C 

R+ = W ■ r SB + W ■ r BC + W ■ r CS + W ' r SB + W ■ r CS m , 

and the average correlation coefficient R in the negative correlation periods: 

w M ■ rf r + w M ■ r" + w c ■ r c + w c ■ r c Rr 

R = M ^ r ^ — . (64) 

W M + W M + W C + W C 

where the weights \sf are the durations of Paleozoic, Mesozoic and Cenozoic respectively: 

w p = 542.0 -251.0 = 291.0 Myr (65) 
w M = 251.0-65.5 = 185.5 Myr (66) 
w c = 65.5 Myr. (67) 

And we denote the difference between R + and R as 

AR = R + -R. (68) 
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We define the average abstract correlation coefficient Q for the positive as well as the negative 
correlation periods as: 

o = wF + ^ + w c £ «' ■ *SJ + i^d + kiD. m 

p=P,M,C 

and the average abstract correlation coefficient Q for the mixtures of positive and negative 
correlation periods as: 

1 

M + W MC 

p=PMC,PM,MC 



Q! 



W PMC + ^PM + W MC 



J] ^-(1^1 + 1^1 + 1^1), 



where the remaining weights w p are: 



(70) 



w 



PMC 



W 



PM 



W 



MC 



542.0 Myr 

542.0 - 65.5 = 476.5 Myr 
251.0 Myr. 



(71) 
(72) 
(73) 



And we denote the difference between Q and Q' as 



AQ = Q-Q'. 



(74) 



We found that the abstract correlation coefficients k™ c |, |r™| and |r^ c | in the mixtures of 



■ fir 



positive and negative periods p = PMC, PM, MC are obviously less than the abstract values |r£ 



flVI 



\rff v \ and \r c \ in the positive or negative periods, namely in the Paleozoic, Mesozoic and Cenozoic 



eras. Therefore, the three climate patterns naturally correspond to the Paleozoic, Mesozoic and 
Cenozoic eras respectively. Based on the data of the first group (n=l) of curves Curve JCC, 
Curve JS L and Curve -BD, we have: 



R + = R + (l 
R- = R~(l 
AR = AR(1 

Q = Qd 
Q' = QV 

AQ = Ag(l 



> (tend to be equal to 1) 

< (tend to be equal to - 1) 

» 

~ 1 (tend to be equal to 1) 

~ (tend to be equal to 0) 

> 



(75) 
(76) 
(77) 
(78) 
(79) 
(80) 



15 



which furthermore shows that the division of three climate patterns CP I, CP II and CP III is 
essential property of the evolutionary earth's spheres. 

Note: These relations are still valid for the other groups of curves (n = 2, 3, 10). 

4.6 The P-Tr extinction was caused by the climate phase reverse between 
CP I and CP II 

We summarize the reasons to explain the P-Tr extinction by the climate phase reverse event as 
follows. 

• Successful explanation of the Sepkoski curve by the tectono-genomic curve based on the 
climate phase reverse event (Fig 1) 

• The climate phase reverse event between CP I and CP II happened at P-Tr boundary (Fig 2a) 

• The sharp peaks of dJCC at the Guadalupian-Lopingian boundary and at the P-Tr boundary 
(Fig 2b) 

• Abnormal climate trend in the Lopingian epoch 

• Different animal extinction patterns at the Guadalupian-Lopingian boundary and at the P-Tr 
boundary. 

4.7 The tectonic curve and the tectonic contribution to the biodiversity 
variation 

The phase of Curve JS L is about the same with the phase of Curve-BD in the Phanerozoic eon. 
And the phase of Curve JCC is about the same with the phase of Curve JiD in the Paleozoic era 
(CP I), while it is about the opposite in the Mesozoic era (CP II) and in the Cenozoic era (CP III). 
Accordingly, we define the associate tectonic curve CurveS ectonicS) by combining the consensus 
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sea level curve and the consensus climate curve as follow (Fig Sle): 

CurveJTectonic J) = [(Curve JS L([P]) + Curve -CC([P])) /2, 

(CurveS L([MC]) - Curve_CC([MC]))/2]. ( } 

We define the tectonic curve Curve -Tectonic with the same standard deviation of the net variation 
biodiversity curve BD: 

Curve _Tectonic = (Curve _TectonicS) - mean(Curve_TectonicS))) ■ a st( i, (82) 

where 

std(BD) 

4 std(Curve -Tectonic J) - mean(CurveS ectonicSf)) 

The tectonic curve Curve -Tectonic represents the tectonic (sea level and climate) contribution 
to the biodiversity evolution. We can calculate the correlation coefficient between the tectonic 
curve and the biodiversity curve in the Paleozoic era or in the Mesozoic and Cenozoic eras: 

r p B = corrcoef(Curve-Tectonic([P]),Curve-BD([P])) 

= 0.421, (84) 

rjjf_ c = corrcoe f '(Curve Sectonic([MC\), Curve _BD([MC])) 

= 0.878. ( } 

Accordingly, we found that the tectonic curve Curve -Tectonic is positively correlated with the 

biodiversity curve Curve SD either in the Paleozoic era or in the Mesozoic and Cenozoic eras. 

5 The genomic contribution to the biodiversity evolution 
5.1 Motivation 

We construct the genomic curve based on the observation of equality between the growth rate kcs 
in genome size evolution and the growth rate k B o in biodiversity evolution. 
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5.2 The overall trend of genome size evolution 



5.2.1 The log-normal distribution of genome size 

We found that the genome sizes of species in a taxon are log-normally distributed in general, which 
were verified in the following 7 taxa (Fig. S2a): 

log(G(A, sp(A))) are normally distributed, (86) 

where G{A, sp(A)) are the genome sizes of all the species sp(A) (sp(A) = 1,2, s m (A)) in the taxon 
A in the genome size databases, and 



A 


= 1 


Diploblostica 


A 


= 2 


Protostomia 


A 


= 3 


Deuterostomia 


A 


= 4 


Bryophyte 


A 


= 5 


Pteridophyte 


A 


= 6 


Gymno sperm 


A 


= 7 


Angio sperm. 



Due to the additivity of normal distribution, the genome sizes of animals, plants, or eukaryotes are 
also log-normal distributed. We obtain the means of logarithm of genome sizes and the standard 
deviations of logarithm of genome sizes as follows: 

GLnjo g W = mean(log(G(A, sp(A)))), (88) 

and 

Gs dJO( ,W = std(log(G(A, sp(A)))\ (89) 

where sp(A) = 1,2, s m (A). Denote G* as the mean logarithm of genome sizes of all the 
contemporary eukaryotes: 

G* = mean(log(G(sp))), (90) 
where sp is all the contemporary eukaryotes in the genome size databases. 

Note: The log-normal distribution of genome size can be demonstrated by the common 
intersection point Q. for the following regression lines (Fig 3b): 

regression line of G meanJog (A') on G sp {A') (91) 
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regression 


line 


of 


GmeanjogiA') + X " GgdJogW) On G sp (A') 


(92) 


regression 


line 


of 


GmeanJogi^') ± X' " G s dj og (A') On G sp (X) 


(93) 


regression 


line 


of 


max(G(A', sp(A'))) on G sp {A') 


(94) 


regression 


line 


of 


min(G(A', sp(A'))) on G sp (A') 


(95) 


regression 


line 


of 


GmeanJogW 011 G P p (A) 


(96) 


regression 


line 


of 


G p , (A) + y ■ G p r , , (A) on G p (A) 

VJ mean Jo g v 7 ~ A sdJog y ' u sp\' K > 


(97) 


regression 


line 


of 


Gl anJog (A) ± X ' ■ G P scUog {A) on G p p {A) 


(98) 


regression 


line 


of 


max(G(A, sp(A))) on G P p (A) 


(99) 


regression 


line 


of 


min(G(A, sp(A))) on G P p (A) 


(100) 



where A = 1,2, 7 for the above 7 taxa, A' = 1,2, 19+53 for 19 animal taxa and 53 angiosperm 
taxa,^ = 1.5677 and^ = 3.1867. The values of G sd j og tend to decline with respect to G sp that is 
proportional to the origin time of taxa (Fig S2b). 

5.2.2 The exponential overall trend of genome size evolution 

We assume the approximate origin times T{A) for the taxa A = 1,2, 7 as follows: 

7X1) = 560.0 Ma 

7X2) = 542.0 Ma, PreCm-Cm 

7X3) = 525.0 Ma 

7/(4) = 488.3 Ma, Cm-0 (101) 

7X5) = 416.0 Ma, S-D 

7X6) = 359.2 Ma, D-C 

7X7) = 145.5 Ma, J-K 

We observed a rough proportional relationship between G p nean log (A) and T(A). Because G p nmn i og (A) 
is the mean genome size of the "contemporary species", we should introduce a new notion (the 
specific genome size) to indicate the mean genome sizes of the "ancient species" in taxa A = 
1,2, 7 at its origin time T(A). Here, we define the specific genome size G p p as: 

G P P (A) = G p eanJog (A) -x ■ G P dJog (A), (102) 

where we let^f = 1.5677 such that the intercept of the regression line of G P p (A) on T(A) is equal to 
G*. We found that G P p (A) is generally proportional to T(A) (Fig. 3a). We define the regression line 
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of G p sp (A) on T(A) as overall trend of genome size curve: 

OT-GS = ko S {-t) + log(N° genome ). (103) 
This equation is equivalent to the exponential overall trend of genome size evolution: 

N genome {t) = N° genome exp(-f/r G5 ), (104) 

where the genome size constant is N2 enome = 2.16 x 10 9 base pairs (bp) and the "e-folding time" 
in genome size evolution is r GS = 256.56 (Myr). The growth rate (namely the slope) of OT-GS is 
k G s = 1/tgs = 0.0038977 Myr 1 . 

Note: The exponential overall trend of genome size evolution obtained in the Phanerozoic eon 
can be extrapolated to the Precambrian period. This extrapolation result according to the value of 
kcs is reasonable to show that the least genome size at 3800 Ma (about the beginning of life) is 
about several hundreds of base pairs (Fig 3d). 



5.3 The agreement between the overall trend of genome size evolution and 
the overall trend of biodiversity evolution 

We found the closely relationship between the genome size evolution and the biodiversity evolution 
(Fig 3d). Both the overall trend of genome size evolution and the overall trend of biodiversity 
evolution are exponential; and the exponential growth rate in the genome size evolution (k GS = 
0.0038977 Myr 1 ) (Fig 3a, 3d) is approximately equal to the exponential growth rate in the 
biodiversity evolution (k BD = 0.0038598 Myr" 1 ) (Fig 2d, 3d): 

k GS * k BD , (105) 

which is equivalent to that the e-folding time in the genome size evolution (r GS = 256.56 Myr) is 
approximately equal to the e-folding time in the biodiversity evolution (tbd = 259.08 Myr): 

t G s ~ r BD - (106) 
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5.4 Explanation of the declining Phanerozoic background extinction rates 



Let ratejjri and rate_ext denote the Phanerozoic biodiversity origination rate and extinction rate 
respectively: 

ratejori : ref. [2], (107) 

rate_ext : ref. [2], (108) 

which agree with each other in general. The difference and the average of them are as follows 
respectively: 

rate - e = (ratejori - rate^exi)/!, (109) 

rate 0+e = (ratejori + rate^ext)/2, (110) 

where rate a - e should agree with d_BD according to their definitions, and rate a+e represents the 
variation of biodiversity in the Phanerozoic eon. The outline of rate a+e indicates the declining 
Phanerozoic background extinction rates [I341 [1351 fl36l PTTl Il38ll . 

We define an essential biodiversity background variation rate by: 

rate ^essential = [amp(\) ■ rate 0+e (l),amp(2) ■ rate 0+e (2), ...,amp(542\) ■ rate 0+e (5421)], (111) 

where 

amp = exp(k GS ■ (-t + 542.0)). (112) 

The outline of ratejessential is generally horizontal (NOT declining). Especially, the peaks of the 
curve rate -essential at P-Tr boundary and at K-Pg boundary are very high, which naturally divide 
the Phanerozoic eon into three climate phases (Fig 2c). 

In the split scenario of biodiversification, we can explain the "declining" background extinction 
rates in the Phanerozoic eon. Firstly, there does not exist a tendency in the essential biodiversity 
background rate curve ratejessential. This essential rate was caused by the random tectonic 
contribution (no tendency) to the biodiversity evolution: 

variation of biodiversity 

ratejessential = . (113) 

tectonic contribution to biodiversity 
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Then, the declining tendency in the observed background extinction or origination rates was caused 
by the genomic contribution to the biodiversity evolution: 

variation of biodiversity 



tectonic contribution+genomic contribution to biodiversity 
It follows that (Fig 2c): 



(114) 



rate 0+e = exp(-k GS ■ (— f + 542.0)) • rate_essential, (115) 

where rate Q+e is declining due to the factor exp(-k GS • {—t + 542.0)). 

The genomic contribution to the biodiversity plays a significant role in the robustness of 
biodiversity evolution: the random tectonic contribution can hardly wipe out all the life on the 
earth thanks for the exponential growth genomic contribution to the biodiversity evolution. 

5.5 Calculating the origin time of taxa based on the overall trend of genome 
size evolution 

5.5.1 The three-stage pattern in Metazoan origination 

We can calculate the origin time of animal taxa according to the linear relationship between the 
origin time and the specific genome size. We obtained the specific genome sizes of the 19 taxa in 
the Animal Genome Size Database (Nematodes, Chordates, Sponges, Ctenophores, Tardigrades, 
Miscellaneous Inverts, Arthropod, Annelid, Myriapods, Flatworms, Rotifers, Cnidarians, Fish, 
Echinoderm, Molluscs, Bird, Reptile, Amphibian, Mammal): 

G animal / i animal \ /^animal / t animal \ , /^animal / i animal \ t\\a\ 
sp \A ) = G meanJog( A ) ~ X ' G sdJog( A )> ( 116 ) 

where \ am,ml = 1,2, 19. We can obtain the origin order of these 19 taxa by comparing their 
specific genome sizes. Hence, we can classify these 19 taxa into Basal metazoa, Protostomia and 
Deuterostomia according to cluster analysis of their specific genome sizes (Data_3). Our result 
supports the three-stage pattern in Metazoan origination based on fossil records fl39l [|40l pT| [|42] 

ma m ea. 
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5.5.2 On angiosperm origination 



Similarly, we can calculate the origin time of angiosperm taxa according to the linear relationship 
between the origin time and the specific genome size. We obtained the specific genome sizes of 
the 53 taxa of angiosperms in the Plant DNA C-value Database (we chose the taxa whose number 
of species is greater than 20 in the calculations): 

/-'angiosperm/ -\angiosperm\ /-.angiosperm / -,angiosperm\ s~iangio sperm/ -,angiosperm\ /i i n\ 

^sp \ A J-*- 7 mean Jog ^ A > X U sdjog ^ A '> \ l 1 ' ) 

where x an8WSperm = 1,2, ...,53. We can obtain the origin order of these 53 taxa by comparing 
their specific genome sizes. Hence, we can classify these 53 taxa into Dicotyledoneae and 
Monocotyledoneae (Data_3). 

Note: The validity of our theory on genome size evolution is supported by its reasonable 
explanation of metazoan origination and angiosperm origination. 

Notation: We denote the mean logarithm genome size, the standard deviation genome size and 
the specific genome sizes by concatenations for all the 19 animal taxa and the 53 plant taxa: 

f~< r (-.animal ^angiosperm-, 

UmeanJog ~ V^meanJoQ' ^ mean lnv J V 11C V 



' mean Jog ~ L ^ mean Jog' ^ mean 

i r /-.animal *~iangios\ 

T sdJog - L ^sdJog' ^sdJog 

< r fanimal fan, 

'sp — L ^ sp ■> S p 



Gr /-'animal /-.angiosperm n /i 1m 

sdjog ~ L TsdJog' U sd loo 

G sp = [ G™ imal , GZ giosperm ]• (120) 



5.6 The phylogenetic tree based on the correlation among genome size 
distributions 



We found that the phylogenetic tree for taxa can be easily obtained based on the correlation 
coefficients among their genome size distributions. We denote the genome size distribution for 
a taxon A by: 

D gs (A, :) = [D gs (A, 1), D gs (A, 2), D gs (A, k), D gs (A, cutoff gs )\ (121) 

where there are D gs (A,k) species in taxon A whose genome size is between (k - 1) • step gs and 
k-step gs , the genome size step step gs = 0.01 picogram (pg) and the genome size cutoff is cutoff gs = 
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2000. Hence, we define the genome size distribution distance matrix M gs (A\, A 2 ) among taxa by: 

M gs (A u A 2 ) = 1 - corrcoef(D gs (A u 0, D gs (A 2 , :)), (122) 
by which, we can draw the phylogenetic tree of the taxa. 

We can obtain the genome size distributions D gs (A, :) and consequently obtain the genome size 
distribution distance matrix M P ^(Ai,A 2 ) among the above 7 taxa as follows: 

M p gs {A x ,A 2 ) = 1 - corrcoef{D p gs {A u 0, D p gs (A 2 , :)), (123) 

where A\,A 2 = 1,2, 7. Hence, we can draw the phylogenetic tree of the 7 taxa based on M gs 
(Fig S2c). 

We can obtain the genome size distributions D a g "" na, (A, :) and consequently obtain the genome 
size distribution distance matrix M a g " imal (A\, A 2 ) among the above 19 animal taxa as follows: 

M animal, ^animal -\animal\ i „„„„„„„ f i r\animai / -tanimal ,\ r^animal/ t animal ,\\ f\ r iA\ 
gs > A 2 ) = 1 _ corrcoe}{D gs (A x , -),O gs (A 2 ,:)), (124) 

where X a ™ mal , A a ™ mal = 1,2, 19. Hence, we can draw the phylogenetic tree of the 19 taxa based 
on M~ l (Fig 3c). 

We can obtain the genome size distributions D a g ® osperm (A, :) and consequently obtain the 
genome size distribution distance matrix Mg" 8losperm (Ai,A 2 ) among the 25 angiosperm taxa (we 
chose 25 angiosperm taxa whose number of species is greater than 50 in the Plant DNA C-value 
database in order to obtain nontrivial distributions) as follows: 

^angiosperm ^angiosperm ^angiosperm^ _ Y—(;Qff(;Qgj'^[) an S' os P erm ^ / \ an S' os P erm jymgiosperm ^angiosperm 

(125) 

where X l ^ glospe}m , \ a ^ los P erm = 1,2, 25. Hence, we can draw the phylogenetic tree of the 25 taxa 
based onM7 ,0 ^™(Fig S2d). 

These phylogenetic trees based on genome size distribution distance matrices generally agree 
with the traditional phylogenetic trees respectively, which is an evidence to show the close 
relationship between the genome evolution and the biodiversity evolution. 



Software: PHYLIP to draw the phylogenetic trees (Neighbor- Joining) in this paper [|46 
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5.7 The varying velocity of molecular clock among taxa 



The growth rates k GS (A) of overall genome size evolution OT taxa {A) for taxa A are not constant, 
though we have an average growth rate k GS for OT-GS . We have an approximate relationship that 
the earlier the origin time T ori (A) is, the slower the growth rate k GS (A) is: 

(k GS (A)-k GS )-T ori (A) = G, (126) 

where the constant G is the difference between the intercept of the overall trend of mean logarithm 
genome size OT meanJog and the intercept of OT-GS . 

5.8 The genomic curve and the genomic contribution to the biodiversity 
evolution 

We define the genomic curve by a straight line with slope k GS and the undetermined intercept b t0 4 ay : 

Curve -Genomic = k GS ■ (-/) + b to d ay , (127) 
which represents the exponential contribution to the biodiversity evolution. 

6 Construction of the tectono-genomic curve 
6.1 The synthesis scheme for the tectono-genomic curve 

The above undetermined intercept of the genomic curve can be defined as: 

btoday = Curve J> epkoski(today) - Curve _Tectonic(today) (128) 
such that Curve _TectonoGenomic(5421) = Curve _S e pkoski{5A2\) . 

We define the tectono-genomic curve by synthesizing the tectonic curve CurveJTectonic and 
the genomic curve Curve -Genomic (Fig 1): 

Curve _T ectonoGenomic = exp(Curve -Tectonic + Curve JJenomic), (129) 
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which agrees very well with the Phanerozoic biodiversity curve Curve _S epko ski: 



Curve _T ectonoGenomic « Curve _S epko ski. 



(130) 



Thus, the Sepkoski curve based on fossil records can be explained by the tectono-genomic curve 
based on climatic, eustatic and genomic data. 

6.2 The driving forces of biodiversity evolution at the molecular level and at 
the species level 

Thus, we have explained the Sepkoski curve in the split scenario. The exponential growth part in 
the Phanerozoic biodiversity evolution was driven by the genome size evolution on one hand, and 
the variation of the the Phanerozoic biodiversity evolution was caused by the Phanerozoic sea level 
fluctuation and climate change on the other hand. 

The successful explanation of the Phanerozoic biodiversity curve CurveJS epko ski shows that 
the driving force of the biodiversity evolution is the tectono-genomic driving force. There are two 
independent tectonic and genomic driving forces in the biodiversity evolution. The first driving 
force originated from the plate tectonics movement at the species level; while the second driving 
force originated from the genome evolution at the molecular level. 

7 The error analysis and reasonability analysis 

7.1 The agreement between the Sepkoski curve and the tectono-genomic 
curve 

7.1.1 The error analysis of the consensus climate curve 

We obtain the first weighted average climate curve C w \ by choosing the corresponding AR(n), 
n = 2, 3, 4 as the weights wl for C 1 , C 2 and C 3 as follows: 



w\ = [AR(2),AR(3),AR(4)]/(AR(2) + AR(3) + AR(4)) 
= [0.3454, 0.1611, 0.4935], 



(131) 
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hence, 



C wl = nondim(w\(\) ■ C 1 + wl(2) • C 2 + wl(3) • C 3 ). 



(132) 



We obtain the second weighted average climate curve C w2 by choosing the corresponding 
correlation coefficients as the weights w2 for C 1 , C 2 and C 3 as follows: 

w2 = [corrcoef(Curve_CC,C l ),corrcoef(Curve_CC,C 2 ), 
corrcoef (Curve _CC, C 3 )]/ 

(corrcoef (Curve J2C, C 1 ) + corrcoef (Curve J2C, C 2 )+ (133) 
+corrcoef '(Curve J2C, C 3 )) 
= [0.4865, 0.2796, 0.2339], 



hence, 



C w2 = nondim(w2(\) ■ C 1 + w2(2) • C 2 + w2(3) • C 3 ). (134) 



We can obtain a weighted average climate curve C w by choosing the average of wl and w2 as 
the weights w for C 1 , C 2 and C 3 as follows: 

w = (wl+w2)/2 

= [0.4159, 0.2204, 0.3637], U ; 

hence, 

C w = nondim(w(\) ■ C 1 + w(2) • C 2 + w(3) • C 3 ), (136) 
which agrees with Curve JCC. 

The weights wl or w2 can be referred to as credibilities for the independent curves C 1 , C 2 
and C 3 . Both of C w \ and C w2 are reasonable estimations of the Phanerozoic climate. So, we can 
consider the zone between C wi and C w2 as the error range of Curve _CC, whose upper range C upper 
and lower range Ci ower are about as follows (Fig Sib): 

Cupper = m UX(C w l,C w 2), (137) 

Ci 0W er = min(C w \,C w2 ). (138) 
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7.1.2 The error analysis of the consensus sea level curve 



We obtain the weighted average sea level curve S w by choosing the corresponding AR(n), n = 

10, 1 1 as the weights w' for S 1 and C 2 as follows: 

W = [AK(10),M(11)]/(AK(10) + AK(11)) 

= [0.4872, 0.5128], 1 ; 

hence, 

S w = nondim(w\Y) • S 1 + w'(2) ■ S 2 ), (140) 

which agrees with Curve SL. 

We can consider the zone between S 1 and S 2 as the error range of CurveSL, whose upper 
range S upper and lower range S !ower are about as follows (Fig Sic): 

Supper = max(S l ,S 2 ), (141) 

Slower = min(S l ,S 2 ). (142) 



7.1.3 The error analysis of the Sepkoski curve 



We can consider the zone between Curve _S _AUGenera and Curve _5 JN ellResolvedGenera as the 
error range of Curve _S epkoski (Fig 1): 

Curve _S MIGenera : ref. [3], (143) 

CurveS JNellResolvedGenera : ref. [3], (144) 

where Curve S JdlGenera is the Phanerozoic biodiversity curve based on all the genera in 
Sepkoski's data and Curve _S -WellResolvedGenera is the Phanerozoic biodiversity curve based 
on well resolved genera in Sepkoski's data. 



7.1.4 The error analysis of the tectono-genomic curve 



In consideration of the error ranges of Curve JZC and CurveSL as well as their phase 
relationships, we define the associate upper tectono-genomic curve Curve _TG_upper_0 and the 
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associate lower tectono-genomic curve Curve JTGJowerJd as follow: 



Curve-TG-upperS) = [(S upper ([P]) + C upper ([P]))/2, (S upper ([MC]) - C lower ([MC]))/2], (145) 

CurveSGJowerD = [(S lowe r([P]) + C lower ([P]))/2, (S !ower ([MC]) - C upper ([MC]))/2] . (146) 

Furthermore, in the similar process and with the same parameters in construction of the 
tectono-genomic curve, we can obtain the upper range and the lower range of the tectono-genomic 
curve as follows (Fig 1): 

Curve_TGMpper = exp(Curve-Genomic+ 

+a st d • (Curve -TG Mp per S) - mean(CurveJTG-upperJ)))), 

Curve-TGJower = exp(Curve-Genomic+ (148) 
+a st d ■ (Curve JTG dower J) - mean(Curve JTGJower _0))). 

7.2 The reasonability of the principal conjectures 

7.2.1 Reasonability of the climate phase reverse based on ^ v (n) 

We can obtain the following 10 groups of curves to describe the Phanerozoic climate, sea level and 
biodiversity: 



n 


= 1 


Curve SL , 


CurveJiD , 


Curve J2C 


n 


= 2 


Curve SL , 


Curve-BD , 


C 1 


n 


= 3 


Curve JiL , 


Curve-BD , 


C- 


n 


= 4 


Curve , 


Curve-BD , 


c 3 


n 


= 5 


Curve SL , 


Curve-BD , 


C w l 


n 


= 6 


Curve SL , 


Curve-BD , 


C W 2 


n 


= 7 


Curve SL , 


Curve-BD , 


c 


n 


= 8 


s l 


Curve-BD , 


CurveJCC 


n 


= 9 


s 2 


Curve-BD , 


CurveJCC 


n 


= 10 


s w , 


Curve_BD , 


Curve-CC 



(149) 



And we can obtain the correlation coefficients r^ v (n) among these groups of curves (Data_2), where 

jj,,v = S , B,C,C ,C ,C , Cwi, C w 2, C w , S , S ,S W (150) 

for the curves CurvSL, Curve -BD, CurveJCC, C 1 , C 2 , C 3 , C wl , C w2 , C w , S 1 , S 2 and S w 
respectively. 
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We can define the corresponding average correlation coefficients for all the 10 groups of curves 
(n = 1, 2, 10) as follows: 

R + (n),R-(n), AR(n), Q(n), Q\n), AQ(n). (151) 

The conclusions on the climate phases CP I, CP II and CP III based on the first group of curves 
(n = 1) still hold for the cases of the other groups of curves (n = 2, 3, 10). Namely, the following 
equations holds in general: 





> 





(152) 


r p BC (n) 


> 





(153) 


r p cs (n) 


> 





(154) 



for CP I, 





> 





(155) 




< 





(156) 




< 





(157) 



for CP II, and 





< 





(158) 


r c BC {n) 


< 





(159) 


r c cs (n) 


> 





(160) 



for CP III. 

Furthermore, we have 

R + (n) > (tend to be equal to 1) (161) 

R'(n) < (tend to be equal to - 1) (162) 

AR(n) » (163) 

Q(n) ~ 1 (tend to be equal to 1) (164) 
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Q'in) ~ (tend to be equal to 0) 
AQ(n) > 



(165) 
(166) 



which shows that the division of three climate phases is an essential property in the evolution rather 
than just random phenomenon in math games. 

The explanation of the P-Tr extinction based on the phase reverse at P-Tr boundary is therefore 
valid regardless the disagreement in the raw data of the Phanerozoic climate and sea level. 
Especially, A7?(l) and A<2(1) are relatively the maximum among these 10 groups of curves, hence 
we chose the optimal first group of curves to describe the Phanerozoic climate, sea level and 
biodiversity throughout this paper. 

The climate system was not stationary when coupling with the other earth's spheres around 
P-Tr boundary. We calculate the correlation coefficients r£ v , where p = P\L,L, L.M.Tr, M\L.M.Tr 
in detail around P-Tr boundary. The curve Curve JCC varies instead in the opposite phase with 
Curve _SL and Curve -BD in Lopingian yet; and it varies instead in the same phase with Curve SL 
and Curve JiD in Lower and Middle Triassic. 

7.2.2 Reasonability of the split scenario 

We summarize the reasons to propose the split scenario in observing the biodiversity evolution as 
follows. 

(1) Evidences to support the close relationship between the genome evolution and the 
biodiversity evolution: 

• Exponential growth in both the genome size evolution and the biodiversity evolution 

• Agreement between genome size growth rate k GS and biodiversity growth rate k BD , namely 

TGS ~ T BD 

• Favorable phylogenetic trees based on M p gs , M"™ mal , M™ giosperm , Mf, M e f 
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• Verification of the three-stage pattern in Metazoan origination and the classification of 
dicotyledoneae and monocotyledoneae in angiosperm origination based on the overall trend 
in genome size evolution 

• Reasonable extrapolation of the overall trend in genome size evolution obtained in 
Phanerozoic eon to the Precambrian period 

• The relationship between phylogenetic trees of species by M ci and the evolutionary tree of 
codons by M CO( i on based on the same matrix A. 

(2) Successful applications of the split scenario: 

• Explanation of the Sepkoski curve by the tectono-genomic curve in the split scenario 

• Error analysis agreement between Curve JS epkoski and Curve _T ectonoGenomic 

• Explanation of the declining Phanerozoic background extinction rates 

• Explanation of the robustness of biosphere in the tremendously changing environment. 



8 The genetic code evolution as the initial driving force in the 
biodiversity evolution 

8.1 The evolutionary relationship between the tree of life and the tree of 
codon 

8.1.1 The codon interval distribution D ci 

We can obtain both the phylogenetic tree of species and the evolutionary tree of 64 codons based 
on the codon interval distributions in the whole genomes. For a certain species a and a certain 
codon n c (n c = 1,2, 64 for 64 codons), we define the "codon interval" I(n c , a, p) as the distance 
between a pair (p) of neighboring codon n c 's in the whole genome sequence. We define the codon 
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interval distribution 



D ci (n c , a, :) = [D ci (n c , a, 1), D ci (n c , a, 2), D ci (n c , a, cutoff ci )]. 



(167) 



as the distribution of all the codon intervals I(n c , a, p) in the whole genome sequence (reading in 
only one direction), where there are D ci {i) pairs of codon n c 's with the distance i (the cutoff of 
distance in the calculations is set as cutoff ci = 1000 bases). For a group of N species, there are 
64 x /V cutoffa-dim vectors D ci {n c , a,:). 

Example: The "GGC" codon interval distribution of the following "genome ao" is 
D ci ("GGC", a , = [0, 0, 1, 3, 5, 1, 0, 0, 0, 0], where cutoffo = 10. 

GGCAUGGCUUGGCAUCGGCAGGCAUGGCAGGCGGCAUGGCAGGCUUGGCAGCA 

And the "GCA" codon interval distribution of the same "genome oro" is D a ("GCA", a , :) = 
[0,0,1,1,2,1,1,0,1,1]. 

GGCAUGGCUUGGCAUCGGCAGGCAUGGCAGGCGGCAUGGCAGGCUUGGCAGCA 

Hence, the correlation coefficient between D a ("GGC", a , :) and D a ("GCA", a , :) is 

corrcoef(D d CGGC\ao, :),D ci ("GCA", a , :)) = 0.7235. 

8.1.2 The codon interval correlation matrix A 

The codon interval correlation matrix A(n c , a,/3) for a group of ,/V species is defined as the 64xNxN 
matrix of the correlation coefficients between pairs of vectors D ci (n c , a) and D ci (n c ,f3): 



8.1.3 Calculating the codon interval distance matrix of species M ci according to A 

We can obtain the N x N codon interval distance matrix M ci (a,/3) of the ./V species by averaging 
the 64 x /V x N correlation coefficients with respect to the 64 codons: 



A(n c , a,P) = corrcoef{D ci (n c , a, :), D ci (n c ,fi, :)). 



(168) 




(169) 
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Hence, we can draw the phylogenetic tree of N species based on M ci . 

The method to obtain phylogenetic trees of species based on the codon interval distance 
matrices is valid not only for eukarya but also for bacteria, archaea and virus. The phylogenetic 
trees of species based on the codon interval distance matrices generally agree with the traditional 
phylogenetic trees respectively, which is also an evidence to show the close relationship between 
the genome evolution and the biodiversity evolution. 

8.1.4 Calculating the distance matrix of codons M codon according to A 

We can obtain the 64 x 64 distance matrix of codons M codon by averaging the 64x NxN correlation 
coefficients with respect to the N species: 

M codon (n c , n' c ) = 1 - corrcoef(A(n c , :, :), A(n' c , :, :)). (170) 

Hence, we can draw the evolutionary tree of 64 codons based on M codon . 

The evolutionary tree of codons based on M co d on agrees with the traditional understanding of 
the genetic code evolution. Thus, we can obtain both the phylogenetic tree of species and the 
evolutionary tree of 64 codons based on the same codon interval correlation matrix A. This is an 
evidence to show the close relationship between the genetic code evolution and the biodiversity 
evolution. The principal rules in the biodiversity evolution may concern the primordial molecular 
evolution. 

8.2 The tree of life and the tree of codon (example 1) 

Based on the genomes of 748 bacteria, 55 archaea, 16 eukaryotes and 133 viruses (GeneBank, 
up to 2009), we can obtain the codon interval correlation matrices A" 11 . For the eukaryotes with 
several chromosomes, the codon interval distributions are obtained by averaging the codon interval 
distributions with respect to the chromosomes of the certain species. Consequently, we can obtain 
the reasonable phylogenetic tree of these speces (Fig S3a) and the reasonable tree of 64 codons 
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(Fig 4a) by calculating M c '! l (a,J3) and M a ", (n c , n' c ) from A aU : 



codon v 

1 M 

(a,P)=l- — Y J A all (n c ,a,{3), (171) 

n c =\ 



and 

MtLn&c, K) = 1 - corrcoef{A"\n c , :, :), A fl "«, :, :)). (172) 

8.3 The tree of life and the tree of codon (example 2) 

Based on the genomes of 16 eukaryotes, we can obtain the codon interval correlation matrices 
A euk . Consequently, we can obtain the reasonable phylogenetic tree of these 16 eukaryotes (Fig 4c) 
and the reasonable tree of 64 codons (Fig S3b) by calculating M™ k (a,fi) and M e c uk don (n c , n' c ) from 
A euk . If there are several chromosomes (chr(a) = 1,2, c m (a)) in the genome of eukaryote a, the 
codon interval distributions of the chromosomes of species a are D e c " k (n c , a, chr(a), :). The codon 
interval correlation matrix is: 

A euk (n c , a, chr(a),j3, chrifi)) = corrcoef(D e f(n c , a, chr(a), :), D e f(n c ,/3, chr(J3), :)). (173) 

Consequently, we can calculating the codon interval distance matrix of species: 

i 64 i i Cm ^ c "'^ 

KfW) = 1 - (a)c (B) £ E A^{n c ,a,chr{a),p,chriP))) (174) 
o<+ n _ x c m \a) c m \p) chr{a)=1 chr{fi)=l 

and the distance matrix of codons: 

MZL^c, n' c ) = 1 - corrcoef(A euk (n c , :,:,:,:), A euk (n' c , :,:,:,:)) (175) 

The phylogenetic tree of eukaryotes by this chromosome average method (for M™ k ) generally 
agrees with the tree by the chromosome average method (for M"! 1 ). 
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8.4 Three periods in genetic code evolution 



We arrange the 64 codons in the "codon_aa" order by considering the codon chronology order 
firstly and considering the amino acid chronology order secondly according to the results in [27]: 

codon chronology: 

(\)GGC,GCC, (2)GUC,GAC, (3)GGG,CCC, (4)GGA,UCC, 

(5)GAG,CUC, (6)GGU,ACC, (7)GCG,CGC, (S)GCU,AGC, 

(9)GCA,UGC, (\0)CCG,CGG, (U)CCU,AGG, {U)CCA,UGG, 

(13)UCG,CGA, (U)UCU,AGA, (\5)UCA,UGA, (\6)ACG,CGU, 

(ll)ACU,AGU, (IS)ACA,UGU, (\9)GAU,AUC, (20)GUG,CAC, 

(2l)CUG,CAG, (22)AUG,CAU, (23)GAA,UUC, (2A)GUA,UAC, 

(25)CUA,UAG, (26)GUU,AAC, (21)CUU,AAG, (2%)CAA,UUG, 

(29)AUA,UAU, (30)AUU,AAU, (3\)UUA,UAA, (32)UUU,AAA, 

amino acid chronology: 

(I) G, (2) A, (3) V, (4) D, (5) P, (6) S, (7) E, (8) L, (9) T, (10) R, 

(II) /, (12)2, (13)iV, (14)*T, (15)//, (16)/% (17)C, (18)M, (19)7, (20) W. 

We define the average correlation curves Hurdle curve and Barrier curve as follows 

Hurdle(a) = mean(M co d on (a, :)), 

Barrieria) = mean({M codon (J3,fi') :\B-a\< n barr and \B' - a\ < n barr }), 
where n barr = 8. 

According to the observations of the certain positions of the three terminal codons in the 
evolutionary tree of codons (Fig 4a, S3b) and the certain shapes of the Hurdle curve and the Barrier 
curve (Fig 4b, S3c), we propose three periods in the genetic code evolution: 

(1) initial period, (2) transition period, and (3) fulfillment period, (180) 

which are separated by the three terminal codons and correspond to the origination of three terminal 
codons respectively. We observe that the curve Barrier begins at a level of Barrier ~ 0.4, then 
overcome a "barrier" of level Barrier ~ 0.5, and at last reach a low place of level Barrier ~ 0.3 
(Fig 4b). Between the initial period and the fulfillment period, we can observe some considerably 
higher values in the curves Barrier and Hurdle, which indicates a "barrier" in the middle period 
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(176) 

(177) 

(178) 
(179) 



of the genetic code evolution. The overall trend of the curve barrier is declining. This "barrier" 
in the curve Barrier corresponds to the narrow palace in the middle of the tree of 64 codons based 
on M codon . 

9 A heuristic model on the coupled earth spheres 

9.1 The strategy of biodiversification 

The robustness of biodiversification was ensured by the genomic contributions, without which 
the biodiversity on the earth can hardly survive the tremendous environmental changes. The 
mechanism of genome evolution is independent from the rapid environmental change during mass 
extinctions, which ensures the continuity of the evolution of life: all the phyla survived from the 
Five Big mass extinctions; more families (in ratio) survived from the mass extinctions than genera. 
The mass extinctions had only influenced some non-fatal aspects of the living system (e.g. wipeout 
of some genera or families), whose influence for the vital or more essential aspects of living system 
(e.g. the advancement aspect) was limited. The living system seems to be able to respond freely 
to any possible environmental changes on the earth. The sustainable development of the living 
system in the high risk earth environment was ensured at the molecular level rather than at the 
species level. 

9.2 The tectonic timescale coupling of earth's spheres 

The three patterns CP I, CP II and CP III in the Phanerozoic eon indicate the tectonic timescale 
coupling of earth's spheres. The driving force in the biodiversity evolution should be explained 
in a tectonic timescale dynamical mechanism. Although the P-Tr mass extinction happened 
rapidly within several 10 4 years, its cause should be explained in a broader context at the tectonic 
timescale. Overemphasis of the impacts of occasional events did not quite touch the core of the 
biodiversity evolution. 
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9.3 A triple pendulum model to explain the climate phase reverse event 

The phase relationship among Curve J D, Curve _SL and Curve JCC can be simulated by a triple 
pendulum model (Fig Sid) with the coupling constants k\, k 2 and a varying coupling k 3 (t) = 
(1 - 6arctan(//?o)/Or/2)) • k 3 : 



This model shows that the climate phase reverse can achieve by just varying the coupling k 3 {t) 
from k 3 (l + e) to k 3 (l - e), e <<c 1. 
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Figure 1: Explanation of the Sepkoski curve by a tectono-genomic curve. 

Curve _T ectonoGenomic generally agrees with Curve JS epkoski not only in overall trends 
but also in detailed fluctuations (including some very detailed fluctuation agreement with 
Curve JS JdlGenerd). The error range of the Sepkoski curve is about between Curve S .AllGenera 
and CurveJS -WellResolvedGenera. The error range of the tectono-genomic curve is about 
between Curve JTGjup 'per and Curve JIG dower . 
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Figure 2: The tectonic contribution to the fluctuations in the biodiversity evolution, a The 

consensus climate curve, the consensus sea level curve and the biodiversification curve. There 
are three climate phases CP I, II, III naturally corresponds to Paleozoic, Mesozoic and Cenozoic 
respectively. Curve -BD generally agrees with Curve JS L. Curve-BD only agrees with Curve JCC 
in the Paleozoic era, but varies oppositely with Curve JCC in the Mesozoic and Cenozoic eras in 
general, b Climate, sea level and biodiversification variation rate curves. We can observe a sharp 
upward peak at GLB and a sharp downward peak at PTB on the curve dJCC. c Explanation 
of the declining Phanerozoic background extinction rates. The overall trend of the essential 
biodiversity background variation rate rate .essential is about horizontal, while the overall trend 
of the origination and extinction rate curves ratejori, ratejext and their average decline due to 
the increasing genomic contribution, d The total biodiversity curve Total-BD is equal to its net 
variation BD plus its overall trend OT-BD. Also, the overall trends of accumulation origination 
and extinction biodiversity curves are exponential. 
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Figure 3: The genomic contribution to the overall trend of the biodiversity evolution, a The 

overall trend in the genome size evolution and its applications: (i) Prediction of origin time of taxa 
in Diploblostica, Protostomia and Deuterostomia indicates three-stage pattern in the metazoan 
origination; (ii) Prediction of origin time of angiosperm taxa differ between Dicotyledoneae and 
Monocotyledoneae. b Proof of the log-normal distribution of genome size in taxa by the common 
intersection point Q. c The phylogenetic tree of animal taxa obtained by M^ s . d Agreement 
between the "e-folding" time tbd in biodiversity evolution and the "e-folding" time tgs in genome 
size evolution. Also, reasonable extrapolation of the overall trend of the genome size evolution 
obtained in the Phanerozoic eon to the Precambrian periods. 
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Figure 4: Relationship between the molecular evolution and the biodiversity evolution, a 

The evolutionary tree of codons obtained by M"'' )don , which agrees with the codon chronolgy. The 
codons in (a) initial period, (b) transition period and (c) fulfilment period are in green, blue and 
red respectively, b The codon distance matrix M" o ' don and its averaging curve Barrier. There was 
a midway high "barrier" {Barrier « 0.5) in the genetic code evolution between the initial period 
and the fulfillment period, c The phylogenetic tree of eukaryotes obtained by their codon interval 
distance matrix M e " k . 
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Figure 5: Fig. Sla The climate curves. 
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Figure 6: Fig. Sib The error range of Curve JCC. 
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Figure 7: Fig. Sic The error range of Curve SL. 
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igure 8: Fig. Sid Simulating climate phase reverse by a triple pendulum model. 
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Figure 9: Fig. Sle The tectonic curve. 
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Figure 10: Fig. S2a Log-normal distributions of genome sizes in taxa. 
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Figure 1 1 : Fig. S2b G sdJog tends to decline with respect to G s 
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Figure 14: Fig. S3a The phylogenetic tree based on M\ 
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Table 1 : Correlation coefficients r£ 
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0.905 


-0.831 




0.703 




0.536 


0.923 


0.957 


0.915 














7 


Curve-BD 


BC W 


-0.006 


-0.517 


-0.889 






-0.763 


0.018 


-0.931 


0.842 


-0.661 


0.494 


-0.654 


1.148 


0.545 


0.361 


0.184 




C„. 


C,S 


0.529 


-0.647 


0.951 




0.617 


-0.921 


0.845 


-0.697 
















s 1 


S l B 


0.506 


0.901 


-0.848 


0.687 


0.431 


0.998 


0.928 


0.888 














8 


Curve-BD 


BC 


0.114 


-0.431 


-0.881 


-0.752 


0.196 


-0.901 


0.205 


-0.638 


0.471 


-0.601 


1.072 


0.511 


0.337 


0.174 




Curve-CC 


CS' 


0.420 


-0.586 


0.915 








0.584 


-0.878 


0.436 


-0.709 
















S' 1 


S 2 B 


0.613 


0.894 


-0.776 




0.611 




0.566 


0.039 


0.906 


0.915 














9 


Curve-BD 


BC 


0.114 


-0.431 


-0.881 


-0.752 


0.196 


-0.901 


0.205 


-0.638 


0.521 


-0.607 


1.128 


0.548 


0.341 


0.207 




Curve-CC 


CS 1 


0.513 


-0.627 


0.909 




0.650 


-0.258 


0.030 


-0.724 
















S w 


S W B 


0.594 


0.905 


-0.830 




0.702 




0.538 


0.915 


0.957 


0.915 














10 


Curve-BD 


BC 


0.114 


-0.431 


-0.881 


-0.752 


0.196 


-0.901 


0.205 


-0.638 


0.516 


-0.611 


1.127 


0.545 


0.343 


0.202 




Curve-CC 


cs w 


0.495 


-0.618 


0.949 




0.659 


-0.902 


0.243 


-0.726 















Table 2: Metazoan origination (G sp = G meanJog ~x • G sdJog ,x = 1-5677) 



No. 


Superphylum 


Taxon 


number of records 


Gjneanjog 


GjsdJog 


G sp 


T ol ; (Ma) 


1 


Protostomia 


Nematodes 


66 


-2.394 


0.93204 


-3.8552 


572.89 


2 


Deuterostomia 


Chordates 


5 


-1.8885 


0.91958 


-3.3301 


566.49 


3 


Diploblostica 


Sponges 


7 


-1.0834 


1.3675 


-3.2272 


565.23 


4 


Diploblostica 


Ctenophores 


2 


-0.010305 


1.6417 


-2.584 


557.39 


5 


Protostomia 


Tardigrades 


21 


-1.2168 


0.7276 


-2.3574 


554.63 


6 


Protostomia 


Miscjnverts 


57 


-0.75852 


0.96321 


-2.2686 


553.54 


7 


Protostomia 


Arthropod 


1284 


-0.078413 


1.2116 


-1.9778 


550 


8 


Protostomia 


Annelid 


140 


-0.14875 


0.9258 


-1.6001 


545.39 


9 


Protostomia 


Myriapods 


15 


-0.54874 


0.66478 


-1.5909 


545.28 


10 


Protostomia 


Flatworms 


68 


0.15556 


1.0701 


-1.522 


544.44 


11 


Protostomia 


Rotifers 


9 


-0.51158 


0.55134 


-1.3759 


542.66 


12 


Diploblostica 


Cnidarians 


11 


-0.16888 


0.69379 


-1.2565 


541.2 


13 


Deuterostomia 


Fish 


2045 


0.23067 


0.6559 


-0.7976 


535.6 


14 


Deuterostomia 


Echinoderm 


48 


0.11223 


0.52794 


-0.71542 


534.6 


15 


Protostomia 


Molluscs 


263 


0.5812 


0.5493 


-0.27994 


529.29 


16 


Deuterostomia 


Bird 


474 


0.32019 


0.13788 


0.10403 


524.61 


17 


Deuterostomia 


Reptile 


418 


0.78332 


0.28332 


0.33916 


521.74 


18 


Deuterostomia 


Amphibian 


927 


2.4116 


1.081 


0.71691 


517.13 


19 


Deuterostomia 


Mammal 


657 


1.1837 


0.2401 


0.80727 


516.03 
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Table 3: Metazoan origination (G' sp = G meanJog - x\ • G sdJog ,x\ = 3.1867) 



No. 


Superphylum 


Taxon 


number of records 


Gjneanjog 


GjsdJog 


G sp 


G sp ' On) 


Ton (Ma) 


3 


Diploblostica 


Sponges 


7 


-1.0834 


1.3675 


-3.2272 


-5.4412 


565.23 


1 


Protostomia 


Nematodes 


66 


-2.394 


0.93204 


-3.8552 


-5.3641 


572.89 


4 


Diploblostica 


Ctenophores 


2 


-0.010305 


1.6417 


-2.584 


-5.2419 


557.39 


2 


Deuterostomia 


Chordates 


5 


-1.8885 


0.91958 


-3.3301 


-4.8189 


566.49 


7 


Protostomia 


Arthropod 


1284 


-0.078413 


1.2116 


-1.9778 


-3.9394 


550 


6 


Protostomia 


Miscjnverts 


57 


-0.75852 


0.96321 


-2.2686 


-3.828 


553.54 


5 


Protostomia 


Tardigrades 


21 


-1.2168 


0.7276 


-2.3574 


-3.5354 


554.63 


10 


Protostomia 


Flatworms 


68 


0.15556 


1.0701 


-1.522 


-3.2545 


544.44 


8 


Protostomia 


Annelid 


140 


-0.14875 


0.9258 


-1.6001 


-3.099 


545.39 


9 


Protostomia 


Myriapods 


15 


-0.54874 


0.66478 


-1.5909 


-2.6672 


545.28 


12 


Diploblostica 


Cnidarians 


11 


-0.16888 


0.69379 


-1.2565 


-2.3798 


541.2 


11 


Protostomia 


Rotifers 


9 


-0.51158 


0.55134 


-1.3759 


-2.2685 


542.66 


13 


Deuterostomia 


Fish 


2045 


0.23067 


0.6559 


-0.7976 


-1.8595 


535.6 


14 


Deuterostomia 


Echinoderm 


48 


0.11223 


0.52794 


-0.7154 


-1.5702 


534.6 


15 


Protostomia 


Molluscs 


263 


0.5812 


0.5493 


-0.2799 


-1.1693 


529.29 


18 


Deuterostomia 


Amphibian 


927 


2.4116 


1.081 


0.71691 


-1.0332 


517.13 


17 


Deuterostomia 


Reptile 


418 


0.78332 


0.28332 


0.33916 


-0.1195 


521.74 


16 


Deuterostomia 


Bird 


474 


0.32019 


0.13788 


0.10403 


-0.1192 


524.61 


19 


Deuterostomia 


Mammal 


657 


1.1837 


0.2401 


0.80727 


0.41857 


516.03 
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Table 4: Metazoan origination (G mean j og ) 



No. 


Superphylum 


Taxon 


number of records 


Gjneanjog 


Gjsdjog 


G sp 


Tori (Ma) 


1 


Protostomia 


Nematodes 


66 


-2.394 


0.93204 


-3.8552 


572.89 


2 


Deuterostomia 


Chordates 


5 


-1.8885 


0.91958 


-3.3301 


566.49 


5 


Protostomia 


Tardigrades 


21 


-1.2168 


0.7276 


-2.3574 


554.63 


3 


Diploblostica 


Sponges 


7 


-1.0834 


1.3675 


-3.2272 


565.23 


6 


Protostomia 


Miscjnverts 


57 


-0.75852 


0.96321 


-2.2686 


553.54 


9 


Protostomia 


Myriapods 


15 


-0.54874 


0.66478 


-1.5909 


545.28 


11 


Protostomia 


Rotifers 


9 


-0.51158 


0.55134 


-1.3759 


542.66 


12 


Diploblostica 


Cnidarians 


11 


-0.16888 


0.69379 


-1.2565 


541.2 


8 


Protostomia 


Annelid 


140 


-0.14875 


0.9258 


-1.6001 


545.39 


7 


Protostomia 


Arthropod 


1284 


-0.078413 


1.2116 


-1.9778 


550 


4 


Diploblostica 


Ctenophores 


2 


-0.010305 


1.6417 


-2.584 


557.39 


14 


Deuterostomia 


Echinoderm 


48 


0.11223 


0.52794 


-0.71542 


534.6 


10 


Protostomia 


Flatworms 


68 


0.15556 


1.0701 


-1.522 


544.44 


13 


Deuterostomia 


Fish 


2045 


0.23067 


0.6559 


-0.7976 


535.6 


16 


Deuterostomia 


Bird 


474 


0.32019 


0.13788 


0.10403 


524.61 


15 


Protostomia 


Molluscs 


263 


0.5812 


0.5493 


-0.27994 


529.29 


17 


Deuterostomia 


Reptile 


418 


0.78332 


0.28332 


0.33916 


521.74 


19 


Deuterostomia 


Mammal 


657 


1.1837 


0.2401 


0.80727 


516.03 


18 


Deuterostomia 


Amphibian 


927 


2.4116 


1.081 


0.71691 


517.13 
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Table 5: Angiosperm origination (G sp = G meanJog ~x • G sdJog ,x = 1-5677) 



1NO. 


Superphylum 


Taxon 


Gjneoji-log 


GsdJtog 






i 
1 


Dicotyledoneae 


Lentibulariaceae 


-l.Ujjz 


U.ooj4y 


-Z.4joZ 


1 7*7 OA 

i / / .yo 


Z 


Monocotyledoneae 


Cyperaceae 


n o n i i 
-U.olZl 1 


n a 1 inn 


1 7710 
-L.I 1 jZ 


1 AO Q< 

loy.oj 


J 


Dicotyledoneae 


Cruciferae 


fl AO 1 OO 

-u.oziyz 


U.ooj j 


1 AOAA 

-i.oyoo 


1 Ay O 1 

ioo\y i 


A 

4 


Dicotyledoneae 


Rutaceae 


n oom 
-U.ZZ1Z1 


n o i/i 1 1 
U.yj41 J 


1 AQ^A 
-l.OoJO 


1 AC "7Q 

loo. /o 


J 


Dicotyledoneae 


Oxalidaceae 


n i a~i7A 
u. iy / /4 


1 1 /l A ^ 

1.1443 


1 ^OA/1 

-i.jyo4 


1 A"7 AO 

lo /.oy 





Dicotyledoneae 


Crassulaceae 


-U.ZOJ /o 


n OOOAO 
U.oZZOo 


i ^^^^ 
-i.jjjj 


1 A"7 1 O 

lo/. iy 


1 


Dicotyledoneae 


Rosaceae 


n /in/1AQ 

-U.4U4oo 


n ao^ i 1 
U.OZjl 1 


1 1Q/1"7 
-l.jo4 / 


1 A^ 1 1 
lOJ. 1 1 


o 
O 


Dicotyledoneae 


Boraginaceae 


Pt OnAA/1 

-U.zUoo4 


n Ay no i 
U.OoUol 


-1.2739 


163.76 


n 

y 


Dicotyledoneae 


Labiatae 


n nno i or\K 

-u.uuziyuj 


U.oUooj 


i nm 
-l.Z /UZ 


i ai n i 
lOJ. / 1 


i n 
1U 


Monocotyledoneae 


Juncaceae 


-U.ZjUjZ 


u.ojoyo 


1 70CO 

-i.zzoy 


1 A "3 J 1 
lOJ.Zl 


1 1 
1 1 


Dicotyledoneae 


Vitaceae 


-u.ouyo / 


n ion/10 

u. jyu4y 


i on 
-1 .ZZZ 


1 A1 1 1 

10j. 1 J 


1 7 

1Z 


Dicotyledoneae 


Cucurbitaceae 


ft OA/1 Q"7 
-U.Z045 / 


u.ou / /y 


1 1 1 T7 
-1.Z1 / / 


i at m 
lOJ.U/ 




Dicotyledoneae 


Onagraceae 


U.U4U545 


U. /oUl o 


1 1 QOO 

-1.1 oZZ 


1 AJ A 1 
10Z.O4 


14 


Dicotyledoneae 


Legummos ae 


n Tin ao 


fl QOAQ/1 

U.oooo4 


i n^nA 
-I.UjUO 


i A i n/i 
101.U4 


1 ^ 

1 j 


Dicotyledoneae 


Myrtaceae 


n ncn i 
-U.J /oUl 


n /i o^ 1 1 
U.4Z3 1 1 


1 f\A A ^ 

- 1 .U44j 


1 An OA 

lou.yo 


i A 
10 


Monocotyledoneae 


Bromeliaceae 


-U.joojo 


fl 7OT30 

u.zyzjz 


i mAA 
-l.UZOO 


i An "7 i 
10U. /4 


1 7 
1 / 


Dicotyledoneae 


Polygon ace ae 


u.zuyoj 


fl "7A 1 "7/1 

U. /Ol /4 


n QQ/i n 
-U.y o4j j 


i An 7 i 
IOU.Zj 


18 


Dicotyledoneae 


Euphorbiaceae 


U. /Zoo / 


1 .0796 


n OA^ A 1 
-U.yojol 


160 


1 Q 
IV 


Dicotyledoneae 


Convolvulaceae 


n ^nn^o 

U. JUUJZ 


n goo 
u.yzo 


n q^/17 
-u.yj4j 


1 SA 

i jy.oo 


on 
ZU 


Dicotyledoneae 


Chenopodiaceae 


n n/i Asno 

-u.u4oouy 


U.jjZo 


n OHIO 
-u.yi jiz 


1 ^O 1A 

i jy. jo 


Zl 


Dicotyledoneae 


Plantaginaceae 


n i <no i 
-U. 1 JUZI 


n /i 8/m 
U.4S4ZZ 


-u.yuyjz 


i jy. j i 


00 

zz 


Dicotyledoneae 


Rubiaceae 


n no /i a 1 1 

-U.U044 1 J 


n ^ 1 ^A^ 
U.J 1 jDj 


n sots 
-u.oyzo 


1 ^O 1 1 

i jy. 1 1 


Zj 


Dicotyledoneae 


Caryophyllaceae 


fl T7AQQ 
U.Z /OOJ 


U.OJooy 


-U. / JJo 


1 ^"7 A A 
1 J / .44 


0/1 
Z4 


Dicotyledoneae 


Amaranthaceae 


U. 1 Joj4 


U.3S1 /O 


-u. / j joy 


1 ^"7 A ") 
1 J / .4Z 


O^ 

Zj 


Dicotyledoneae 


Malvaceae 


n io^ 1 1 

u.jyji / 


n a ~i i no 

u.4/ iuy 


n iaha 
-U. j4jjo 


1 ^.7 A 1 
1 jZ.41 


OA 
ZO 


Monocotyledoneae 


Zingiberaceae 


fl 7AQ 1 o 
U.Z4oiy 


U.jOj 1 / 


n n i i ^ 
-U. jZI 1 j 


1 ^7 \ A 

1 jZ. 14 


T7 
Z / 


Monocotyledoneae 


Iridaceae 


1 1/IOQ 

i.j4zy 


i.u4yi 


n in i ts 
-U.jUI /o 


i ji.y 


28 


Dicotyledoneae 


Umbelliferae 


r» -t i nm 
U. / IUUj 


0.6235 


-U.Zo /4Z 


151.49 


OQ 

zy 


Dicotyledoneae 


Solanaceae 


U. /oU34 


U.OOJOJ 


n TAl^T 

-U.zojjz 


1^1 A A 
1 jl .44 


in 
jU 


Monocotyledoneae 


Orchidaceae 


1 /lOAl 

1 .4U0 J 


I.Ujj 1 


n HAHQA 

-U.Z4 /o4 


1 jI.Zj 


J 1 


Monocotyledoneae 


Araceae 


1 ^ 1 7 A 
l.J 1 /4 


inn 
1.U1Z 


n nAO i 
-u.uoy 1 jz 


1 AO 0"7 

i4y.u / 


10 

3Z 


Dicotyledoneae 


Pap ave race ae 


n oionA 

u.yjzuo 


n a 1 on 
U.oiyjZ 


-U.Ujooj4 


1 A Q 7 
140. / 


j J 


Dicotyledoneae 


Compositae 


1 (Y7 1 1 

1 .U/41 


n "7moA 
U. /U /Zo 


n f\i 1 a^7 
-U.Uj4oj / 


1 1 Q A^ 
140.OJ 


1/1 
j4 


Monocotyledoneae 


Gramineae 


1 /inm 
1 .4UUZ 


n Q/I/ITA 

U.o44 /O 


U.U / joy4 


1 A7 1 
14 /. J 


jj 


Dicotyledoneae 


Cactaceae 


n noo 1 1 
U.yool j 


n cm; 1 
U.j /Zjl 


n nonAns 
U.UyUoUo 


1 /I "7 1 ") 
14 /. IZ 


1A 
JO 


Monocotyledoneae 


Palmae 


1 1 TOT 

1 . 1 ZZZ 


n Ai/i oq 
U.dj4oo 


n 1 TAO 1 

u.izoyi 


1 /I A AQ 
140.00 


17 


Dicotyledoneae 


Passifloraceae 


u. jzzuy 


U.ZZ4/Z 


n 1 AOTO 

u. ioy /y 


1 1A 1 ^ 
140. 1 J 


38 


Dicotyledoneae 


Orobanchaceae 


1.12 


0.54393 


n OA"7JA 
U.ZO /Zo 


144.97 


1Q 

jy 


Dicotyledoneae 


Cistaceae 


n cson^ 
u.ooyuj 


U. JU4JO 


n a 1 1 ^ 

U.4 1 1 JJ 


1/1101 
14J.Z 1 


40 


Monocotyledoneae 


Asparagaceae 


2.0053 


0.78802 


0.76991 


138.84 


41 


Dicotyledoneae 


Asteraceae 


1.8795 


0.67031 


0.82863 


138.12 


42 


Dicotyledoneae 


Ranunculaceae 


2.0285 


0.72517 


0.8916 


137.35 


43 


Monocotyledoneae 


Agavaceae 


1.6207 


0.4537 


0.90941 


137.13 


44 


Monocotyledoneae 


Hyacinthaceae 


2.3635 


0.69028 


1.2814 


132.6 


45 


Dicotyledoneae 


Loranthaceae 


2.3797 


0.68478 


1.3062 


132.3 


46 


Monocotyledoneae 


Commelinaceae 


2.5322 


0.64196 


1.5258 


129.62 


47 


Monocotyledoneae 


Amaryllidaceae 


2.9085 


0.5811 


1.9975 


123.86 


48 


Monocotyledoneae 


Xanthorrhoeaceae 


2.7036 


0.4 


2.0765 


122.9 


49 


Monocotyledoneae 


Asphodelaceae 


2.8054 


0.33968 


2.2729 


120.51 


50 


Monocotyledoneae 


Alliaceae 


2.9051 


0.39078 


2.2924 


120.27 


51 


Dicotyledoneae 


Paeoniaceae 


2.957 


0.28164 


2.5155 


117.55 


52 


Monocotyledoneae 


Liliaceae 


3.5678 


0.63278 


2.5757 


116.81 


53 


Monocotyledoneae 


Aloaceae 


2.9724 


0.2203 


2.627 


116.19 
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Table 6: Angiosperm origination (G' sp = G meanJog - x\ • G sdJog ,x\ = 3.1867) 



No. 


Superphylum 


Taxon 


GjneanJog 


GsdJog 




G'sp 


Tori (Ma) 


1 


Dicotyledoneae 


Lentibulariaceae 


-1.0532 


A O 1 A A 

0.88349 


-2.4382 


-3.8686 


177.96 


5 


Dicotyledoneae 


Oxalidaceae 


A 1 f\11 A 

0.19774 


1.1445 


-1.5964 


-3.4494 


167.69 


4 


Dicotyledoneae 


Rutaceae 


-0.22121 


0.93413 


-1.6856 


-3.198 


168.78 


6 


Dicotyledoneae 


Crassulaceae 


-0.265 /o 


0.82268 


-1.5555 


-2.8874 


167.19 


3 


Dicotyledoneae 


Cruciferae 


-0.62192 


0.6855 


-1.6966 


-2.8064 


168.91 


2 


Monocotyledoneae 


Cyperaceae 


-0.8121 1 


0.61307 


-1.7732 


-2.7658 


169.85 


18 


Dicotyledoneae 


hupnorbiaceae 


a i o /; o i 

0. /2687 


1.0796 


-0.96561 


-2.7135 


160 


9 


Dicotyledoneae 


Labiatae 


n (\t\i i (iik 


A OAOO 1 


-1.2702 


-2.5797 


163.71 


1 A 

14 


Dicotyledoneae 


Leguminosae 


A lOOiCO 

o.jjyoo 


("A o o /; o i 


-1.0506 


-2.4864 


161.04 


19 


Dicotyledoneae 


Convolvulaceae 


0.50052 


0.928 


-0.9543 


-2.4567 


159.86 


13 


Dicotyledoneae 


Onagraceae 


r\ t\ a r\ O A O 

0.040848 


A "7 O A 1 O 

0.78018 


-1.1822 


-2.4454 


162.64 


7 


Dicotyledoneae 


Rosace ae 


A A(\A£^<D 

-0.40468 


0.6251 1 


-1.3847 


-2.3967 


165.1 1 


8 


Dicotyledoneae 


Boraginaceae 


A 1 A/l/i A 

-0.20664 


A £L O AO 1 

0.68081 


-1.2739 


-2.3762 


163.76 


10 


Monocotyledoneae 


Juncaceae 


-0.230 3 2 


0.63698 


-1.2289 


-2.2602 


163.21 


17 


Dicotyledoneae 


Polygonaceae 


a ~>nno c 

0.20985 


0.76174 


-0.98433 


-2.2176 


160.23 


12 


Dicotyledoneae 


Cucurbitaceae 


a i c a o i 

-0.2648/ 


A /LAT7A 

0.60779 


-1.2177 


-2.2017 


163.07 


27 


Monocotyledoneae 


Iridaceae 


1.3429 


1.0491 


A ~) A 1 1 

-0.30178 


-2.0003 


151.9 


30 


Monocotyledoneae 


Orchidaceae 


1.4063 


1.0551 


A 1 A 1 A 

-0.24/84 


-1.956 


151.25 


1 1 


Dicotyledoneae 


Vitaceae 


a iJAnQ'7 

-U.6Uyo / 


a i on i n 
0.jV04y 


-1.222 


-1.8542 


163. 13 


zi 


Dicotyledoneae 


Caryophyllaceae 


O.z /6o j 


a /rr o/:n 

U.6!>o6y 


-0.7558 


-1.8222 


157.44 


20 


Dicotyledoneae 


Chenopodiaceae 


a A A £.Q no 

-0.046809 


0.5526 


A A 1 1 1 1 

-0.91312 


-1.8078 


159.36 


15 


Dicotyledoneae 


Myrtaceae 


A T "70 A 1 

-0.3 /801 


0.4251 1 


-1.0445 


-1.7327 


160.96 


22 


Dicotyledoneae 


Rubiaceae 


A AO A A 1 O 

-0.084413 


0.51565 


A O A 1 O 

-0.8928 


-1.7276 


159.1 1 


31 


Monocotyledoneae 


Araceae 


1.5174 


1.012 


-0.069152 


-1.7075 


149.07 


24 


Dicotyledoneae 


Amaranthaceae 


0.15834 


0.58176 


-0. Z5369 


-1.6956 


157.42 


21 


Dicotyledoneae 


Plantaginaceae 


-0.15021 


0.48422 


A AAA~>^ 

-0.90932 


-1.6933 


159.31 


16 


Monocotyledoneae 


Bromeliaceae 


-0.5683s 


0.29232 


-1.0266 


-1.4999 


160.74 


29 


Dicotyledoneae 


Solanaceae 


A 1 0AT /I 

0.78034 


0.66585 


-0.26352 


-1.3415 


151.44 


34 


Monocotyledoneae 


uramineae 


1 .4002 


A O A A 1 C 

0.84476 


A AT cond 

0.0/5894 


-1.2918 


147.3 


28 


Dicotyledoneae 


Umbelliferae 


("1 "7 1 AA"> 

0.71003 


0.6235 


-0.26742 


-1.2769 


151.49 


33 


Dicotyledoneae 


Compositae 


1.0741 


A 1 ATI/; 

0.70726 


a m a 

-0.034657 


-1.1797 


148.65 


25 


Dicotyledoneae 


Malvaceae 


0.39517 


A Al t AA 

0.47109 


-0.34336 


-1.1061 


152.41 


32 


Dicotyledoneae 


Papaveraceae 


A A A/; 

0.93206 


0.61932 


A A "5 O O C A 

-0.038854 


-1.0415 


148.7 


26 


Monocotyledoneae 


Zingiberaceae 


0.24819 


0.36317 


-0.321 15 


-0.9091 


152.14 


36 


Monocotyledoneae 


Palmae 


1.1222 


A H'} AOO 

0.63488 


0.12691 


-0.901 


146.68 


35 


Dicotyledoneae 


Cactaceae 


A A O O 1 1 

0.98813 


0.57251 


A AAA/" AO 

0.090608 


-0.8363 


147.12 


38 


Dicotyledoneae 


Orobanchaceae 


1.12 


0.54393 


0.26726 


-0.6133 


144.97 


40 


Monocotyledoneae 


Asparagaceae 


2.0053 


A "7 0A1 

0.78802 


A "7/1 A A 1 

0.76991 


-0.5059 


138.84 


A 1 
4z 


Dicotyledoneae 


Ranunculaceae 


Z.UzoD 


U. /zo 1 / 


n on i a 


-U.zoz4 


1 J 1 .5 J 


41 


Dicotyledoneae 


Asteraceae 


1.8795 


0.67031 


0.82863 


-0.2566 


138.12 


37 


Dicotyledoneae 


Passifloraceae 


0.52209 


0.22472 


0.16979 


-0.194 


146.15 


39 


Dicotyledoneae 


Cistaceae 


0.88905 


0.30458 


0.41155 


-0.0816 


143.21 


44 


Monocotyledoneae 


Hyacinthaceae 


2.3635 


0.69028 


1.2814 


0.16378 


132.6 


43 


Monocotyledoneae 


Agavaceae 


1.6207 


0.4537 


0.90941 


0.17489 


137.13 


45 


Dicotyledoneae 


Loranthaceae 


2.3797 


0.68478 


1.3062 


0.19751 


132.3 


46 


Monocotyledoneae 


Commelinaceae 


2.5322 


0.64196 


1.5258 


0.48647 


129.62 


47 


Monocotyledoneae 


Amaryllidaceae 


2.9085 


0.5811 


1.9975 


1.05671 


123.86 


48 


Monocotyledoneae 


Xanthorrhoeaceae 


2.7036 


0.4 


2.0765 


1.42892 


122.9 


52 


Monocotyledoneae 


Liliaceae 


3.5678 


0.63278 


2.5757 


1.55132 


116.81 


50 


Monocotyledoneae 


Alliaceae 


2.9051 


0.39078 


2.2924 


1.6598 


120.27 


49 


Monocotyledoneae 


Asphodelaceae 


2.8054 


0.33968 


2.2729 


1.72294 


120.51 


51 


Dicotyledoneae 


Paeoniaceae 


2.957 


0.28164 


2.5155 


2.0595 


117.55 


53 


Monocotyledoneae 


Aloaceae 


2.9724 


0.2203 


2.627 


2.27037 


116.19 
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Table 7: Angiosperm origination (G mean j og ) 



JNO. 


Superphylum 


Taxon 


Gjfiectfi-log 


GsdJtog 




l ori (Ma) 


i 
1 


Dicotyledoneae 


Lentibulariaceae 


-l.UjjZ 


U.ooJ4y 


1 A 1C1 

-Z.4joZ 


1 T*7 (1A 

i / / .yo 


Z 


Monocotyledoneae 


Cyperaceae 


n cm i 

"U.OlZl 1 


fi A 1 7fi7 
U.Ol JU/ 


1 *7*7'20 
-1. / / JZ 


loy.oj 


1 
J 


Dicotyledoneae 


Cruciferae 


fi A 7 1 07 

-u.oziyz 


fi ^o^*; 
U.Ooj j 


-i.oyoo 


los.y i 


i i 
1 1 


Dicotyledoneae 


Vitaceae 


-u.ouys / 


fi IQfi/IQ 

u. jyu4y 


-1 .zzz 


1 Ai 1 1 

lOJ. 1 J 


1 A 
10 


Monocotyledoneae 


Bromeliaceae 


fi ^ACJC 

-U.joojo 


fi 000*20 

u.zyzjz 


-l.UZOO 


1 Afi "7 1 
lOU. /4 


1 


Dicotyledoneae 


Rosaceae 


fi /Ifi/IAQ 
-U.4U400 


fi /io*; i 1 
U.oZjl 1 


-l.Jo4 / 


i At; i i 
lOJ. 1 1 


1 c 

Ij 


Dicotyledoneae 


Myrtaceae 


-U.j /oUl 


n /1 7 *i i i 
U.4Zj 1 1 


1 fi/1 A *J 

- 1 .U44j 


1 Afi OA 

lou.yo 


/: 
O 


Dicotyledoneae 


Crassulaceae 


fi 7A*;7Q 
-U.ZOJ /O 


U.oZZoo 


-1.5555 


167.19 


1 7 
1Z 


Dicotyledoneae 


Cucurbitaceae 


fi 7 A/1 Q"7 

-U.Zo4o / 


fi AfiOOO 

u.ou/ /y 


1 1 1 *7*7 
-1.Z1 / / 


1 fi.1 fi*7 
lOJ.U / 


i fi 

1U 


Monocotyledoneae 


Juncaceae 


fi OQfilO 

-U.ZjUjZ 


u.ojoyo 


i nco 

-i.zzoy 


1 f.1 O 1 
lOJ.Zl 


A 

4 


Dicotyledoneae 


Rutaceae 


fi oom 
-U.ZZ1Z1 


fi O l/I 1 1 

u.yj4i j 


1 ao*;a 
-I.Oojo 


1 AC "70 
loo. la 


o 
o 


Dicotyledoneae 


Boraginaceae 


fi ofiAA/i 
-U.ZU004 


fi Acne i 
U.ooUol 


-i.z / jy 


lOJ. /O 


Zl 


Dicotyledoneae 


Plantaginaceae 


fi 1 *CfiO 1 

-U.l jUZI 


fi /I Q/1 77 

U.4o4ZZ 


fi ofion 

-u.yuyjz 


1 Jy.i 1 


77 

ZZ 


Dicotyledoneae 


Rubiaceae 


fi HQ /I /I 1 1 

-U.Uo441 j 


U.J 1 jOj 


fi CQTO 

-u.syzo 


i jy. 1 1 


7fi 
ZU 


Dicotyledoneae 


Chenopodiaceae 


-u.u4oouy 


fi *;*;ia 
U.jjZo 


fi Q1 1 1 O 

-u.yi jiz 


1 *nQ 1A 

i jy. jo 


n 

y 


Dicotyledoneae 


Labiatae 


-u.uuziyuo 


U.oUooj 


1 0*7fiO 

-i.z /uz 


1 fi.1 7 1 
lOJ. / 1 


1 1 


Dicotyledoneae 


Onagraceae 


fi fi/lfiC 1 Q 


fi 7Qfi1 

U. /oUl o 


-1.1 oZZ 


1 A") A 1 
10Z.04 


24 


Dicotyledoneae 


Amaranthaceae 


0. 15834 


0.58176 


-0.75369 


157.42 


J 


Dicotyledoneae 


Oxalidaceae 


fl 1 0*7*7/1 

u. iy / /4 


1 1 A A *\ 
1 . 144J 


- 1 .jyo4 


1 A7 AQ 

io/. oy 


1 7 
1 / 


Dicotyledoneae 


Polygon ace ae 


fi OfiQQ< 

u.zuyoj 


fi 7 A 1 "7/1 
U. /Ol /4 


fi QO/1 11 

-U.y o4j 3 


1 Afi 71 
lOU.Zj 


7A 
ZO 


Monocotyledoneae 


Zingiberaceae 


fi 7/1 Q 1 O 

U.z4oiy 


fi 7 A7 1 "7 
U. JOJ 1 / 


fi ni i ^ 
-U. jZI 1 j 


1 <0 1/1 
1 jZ. 14 


77 
Zj 


Dicotyledoneae 


Caryophyllaceae 


fl 77AS7 
U.Z /OoJ 


fi A*\GAQ 


fi T^'sS 
-U. / jjo 


1 J / .44 


1 /I 
14 


Dicotyledoneae 


Legummos ae 


fi IIOAO 


fi OC^Cyl 
U.OOO04 


-I.UjUo 


1 A 1 fi/1 
101.U4 


7 *" 

Zj 


Dicotyledoneae 


Malvaceae 


n "jo*; i *7 

u.jyji / 


fi /1*7 1 fiO 

u.4/ luy 


fi 1A11& 

-U. j4jjO 


1^7 ,-1 1 
1 jZ.41 


1 Q 

iy 


Dicotyledoneae 


Convolvulaceae 


fi *ifin*"7 
U. jUUjZ 


fi QOQ 

u.yzo 


fi o*;/i "j 

-u.yj4j 


i jy.JSo 


77 


Dicotyledoneae 


Passifloraceae 


fi COOfiQ 

u. jzzuy 


fi 10/1*70 
U.ZZ4/Z 


fi 1 A 0*70 

u. ioy /y 


1 1A It: 
14o. 1 J 


7Q 

Zo 


Dicotyledoneae 


Umbelliferae 


fi 7 1 nr\i 
U. / 1UUJ 


fi ATI*; 
U.OZjj 


fi 1^7 /lO 
-U.Zo /4Z 


1 *N 1 /IQ 

i ji.4y 


18 


Dicotyledoneae 


Euphorbiaceae 


fi Tl/^OT 

U. /zoo / 


1.0796 


fi n^*;^i 

-u.yojoi 


160 


7Q 

zy 


Dicotyledoneae 


Solanaceae 


fi 7Qfi7/1 

U. /oU34 


U.oojoj 


fi OAKO 

-U.ZojjZ 


i *; i /i /i 
1 jl .44 


7Q 

jy 


Dicotyledoneae 


Cistaceae 


fi QQOfi< 

u.ooyuj 


fi "Jfi/1 

U. jU4jo 


(\ A] 1 f \ f ; 
U.41 1 jj 


1/17 7 1 

14J.Z1 


77 

JZ 


Dicotyledoneae 


Papaveraceae 


n Q77fiA 

u.yjzuo 


fi a 1 on 

u.oiyjz 


fi fiQCQ^-l 

-U.Ujooj4 


1 A O 7 
14o. / 


7 *" 

JJ 


Dicotyledoneae 


Cactaceae 


n noo i 7 

u.yooi j 


fi <*70'* 1 

U.J /Zjl 


fi fiOfiAfiQ 

u.uyuouo 


1/17 17 
14 /. IZ 


7 7 

J J 


Dicotyledoneae 


Compositae 


1 fi"7 1 1 

1 .U/41 


fi *7fi"7TA 

u. /u /zo 


fi fi Q -1 ^.*v*7 

-U.Uj4oj / 


1 1 Q A^ 
148. OJ 


1Q 

JO 


Dicotyledoneae 


Orobanchaceae 


1 1 1 
1. 1Z 


fi *I/1 lOT 

u. j4jyj 


fi OA'70A 

U.Zo /Zo 


1/1/1 Q7 

144. y / 


JO 


Monocotyledoneae 


Palmae 


1. 1ZZZ 


fi £.1 1 OQ 

U.oj4oo 


fi 1 1AO 1 

u.izoy i 


1 A A AQ 
140.UO 


77 

Z / 


Monocotyledoneae 


Indaceae 


1 7/17Q 

i.j4zy 


i.U4y i 


fi Ifil "7Q 

-U.jUI /8 


1 *" 1 O 

i ji.y 


1A 

j4 


Monocotyledoneae 


Gramineae 


1 .4UUZ 


fi C/1/1*7A 

U.o44 /O 


fi fiT^Cfi 1 

u.u / joy4 


1/17 7 
14 /. J 


30 


Monocotyledoneae 


Orchidaceae 


1 .4063 


1.0551 


fi A HQ A 

-U.Z4 /o4 


151.25 


J 1 


Monocotyledoneae 


Araceae 


1 *\ 1 7/1 
1 . J 1 /4 


1 fi 1 T 
1 .U1Z 


fi fiAO 1 *\0 

-u.uoy i jz 


1 /1Q fi7 

i4y.u / 


43 


Monocotyledoneae 


Agavaceae 


1.6207 


0.4537 


0.90941 


137.13 


41 


Dicotyledoneae 


Asteraceae 


1.8795 


0.67031 


0.82863 


138.12 


40 


Monocotyledoneae 


Asparagaceae 


2.0053 


0.78802 


0.76991 


138.84 


42 


Dicotyledoneae 


Ranunculaceae 


2.0285 


0.72517 


0.8916 


137.35 


44 


Monocotyledoneae 


Hyacinthaceae 


2.3635 


0.69028 


1.2814 


132.6 


45 


Dicotyledoneae 


Loranthaceae 


2.3797 


0.68478 


1.3062 


132.3 


46 


Monocotyledoneae 


Commelinaceae 


2.5322 


0.64196 


1.5258 


129.62 


48 


Monocotyledoneae 


Xanthorrhoeaceae 


2.7036 


0.4 


2.0765 


122.9 


49 


Monocotyledoneae 


Asphodelaceae 


2.8054 


0.33968 


2.2729 


120.51 


50 


Monocotyledoneae 


Alliaceae 


2.9051 


0.39078 


2.2924 


120.27 


47 


Monocotyledoneae 


Amaryllidaceae 


2.9085 


0.5811 


1.9975 


123.86 


51 


Dicotyledoneae 


Paeoniaceae 


2.957 


0.28164 


2.5155 


117.55 


53 


Monocotyledoneae 


Aloaceae 


2.9724 


0.2203 


2.627 


116.19 


52 


Monocotyledoneae 


Liliaceae 


3.5678 


0.63278 


2.5757 


116.81 
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Table 8: Origination order 



Superphylum 


T„ri (Ma) 


G mean Jog 


G sdjog 


G sp 




Gin in 


Diploblostica 


560 


-0.4731 


1.095 


-2.1898 


1.1506 


-2.8134 


Protostomia 


542 


-0.10229 


1.2158 


-2.0083 


4.1685 


-3.912 


Deuterostomia 


525 


0.87752 


1.0869 


-0.82636 


4.8891 


-2.8134 


bryophyte 


488.3 


-0.63576 


0.54685 


-1.4931 


2.0757 


-1.772 


pteridophyte 


416 


1.7359 


1.6606 


-0.86744 


4.2861 


-2.4079 


gymnosperm 


359.2 


2.8263 


0.46055 


2.1043 


3.5835 


0.81093 


angiosperm 


145.5 


0.96878 


1.2681 


-1.0193 


5.0252 


-2.8134 


Protist 




-1.5532 


1.6488 


-4.1381 


2.9755 


-7.3475 


Eubacteria 




-5.8238 


0.57889 


-6.7313 


-4.5865 


-8.7269 


Ai'chaea 




-6.02 


0.50451 


-6.811 


-4.7404 


-6.9616 
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